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Oh Abstract. Nongaussian statistics are a powerful discriminant between inflationary models, 

O particularly those with noncanonical kinetic terms. Focusing on theories where the La- 

+J grangian is an arbitrary Lorentz-invariant function of a scalar field and its first derivatives, 

Ci we review and extend the calculation of the observable three-point function. We compute 

the "next-order" slow-roll corrections to the bispectrum in closed form, and obtain quan- 
titative estimates of their magnitude in DBI and power-law /c-inflation. In the DBI case 
\^ our results enable us to estimate corrections from the shape of the potential and the warp 

^ factor: these can be of order several tens of percent. We track the possible sources of large 

■nj- logarithms which can spoil ordinary perturbation theory, and use them to obtain a general 

formula for the scale dependence of the bispectrum. Our result satisfies the next-order version 
of Maldacena's consistency condition and an equivalent consistency condition for the scale 
j—{ dependence. We identify a new bispectrum shape available at next-order, which is similar 

^ to a shape encountered in Galileon models. If /nl is sufficiently large this shape may be 

independently detectable. 
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1 Introduction 



In the last decade, our picture of the early universe has become numerically precise — a 
quantitative revolution made possible by analysis of the cosmic microwave background (CMB) 
radiation. The CMB is usually interpreted as a relic of an earlier hot, dense, primordial era [1] . 
During this era the universe is believed to have been filled with an almost-smooth plasma, 
in which tiny perturbations were seeded by an unknown mechanism. Eventually, photons 
decoupled from the cooling plasma to form the CMB, while matter particles collapsed under 
gravity to generate bound structures. 

The origin of these perturbations is a matter of dispute. A possible candidate is the 
inflationary scenario, in which an era of quasi-de Sitter expansion smoothed the universe 
while allowing quantum fluctuations to grow and become classical on superhorizon scales. In 
its simplest implementation, inflation forecasts an approximately scale-invariant, Gaussian 
distribution of perturbations [2, 3], characterized entirely by their two-point correlations. 
Although falsifiable, these predictions are in agreement with present-day observations [4-6]. 

The increasing sophistication of CMB experiments has suggested that nontrivial three- 
and higher n-point correlations could be detected [7, 8], and in the medium-term future 
it is possible that competitive constraints will emerge from the nongaussian statistics of 
collapsed structures [9]. In principle, valuable information is encoded in each n-point function. 
In practice, extracting information from the four-point function is already computationally 
challenging [10, 11], and it is unclear whether useful constraints can be obtained for the five- 
and higher n-point functions. For this reason attention has focused on the amplitude of 
three-point correlations, /nl [7, 8]. For a Gaussian field these correlations are absent and 
/nl = 0. 

In simple models both the de Sitter expansion and quantum fluctuations originate from 
a single scalar field, (j), although more complicated possibilities exist. ^ But whatever combi- 
nation of ingredients we choose, the perturbations inherit their statistical properties from a 
microphysical Lagrangian by some mapping of the correlation functions of a field such as (j) 
[13-15] . The two-point function for a scalar field with canonical kinetic term was calculated 
during the early days of the inflationary paradigm [2, 16-21] , followed by partial results for the 
three-point function [22-25] . The complete three-point function was eventually computed by 
Maldacena [26], and later extended to multiple fields [27] and four-point correlations [28, 29]. 

Maldacena's calculation showed that /nl would be unobservably small in the simplest 
model of inflation — of order the tensor-to-scalar ratio, r, which is constrained by observation 
to satisfy r < 0.2 [6]. However, it was clear that information about the Lagrangian operators 
responsible for generating fluctuations in more exotic theories could be encoded in the three- 
and higher n-point functions. Which operators could be relevant? Creminelli argued that if 
the dominant kinetic operator for the background field was the canonical term (Vcj))'^, where 
V is the spacetime covariant derivative, then the three-point correlations of its perturbations 
were effectively indistinguishable from Maldacena's simplest model [30]. However, in theories 
where other operators make significant contributions, the amplitude of three-point correla- 

^For a review of the various generation mechanisms which have been proposed to date, see Lyth & Liddle 
[12]. 
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tions can be radically different. This gave valuable diagnostics for more complicated scenarios 
such as "ghost inflation" [31, 32] or models based on the Dirac-Born-Infeld (DBI) action [33]. 
In these examples, and many others, the scalar field action controlling the dynamics of both 
background and perturbations can be written 

S=l- [d^x^\R + 2P{X,cl))\, (1.1) 



2 „ 

where X = —g°'^\7a(l)^b4'i the spacetime Ricci scalar is R, and P{X, (p) is an arbitrary func- 
tion. More recently, diagnostic signatures for related models based on "Galileon" actions 
have been obtained [34-38]. We describe models in which the background kinetic structure 
is dominated by (V(/))^ as canonical, and more general examples of the form (1.1) as non- 
canonical. Eq. (1.1) was suggested by Armendariz-Picon et al. [39], who described it as 
'A;-inflation.' The corresponding two-point function was obtained by Garriga & Mukhanov 
[40]. Three-point correlations were studied by Gruzinov in a decoupling limit where mixing 
with gravity could be neglected [41]. Gravitational interactions were included in Ref. [42]. 
The bispectrum was given in generality by Chen et al. [43], followed by extensions to multiple 
fields [44-46]. The four-point function was computed by several authors [47-52]. 

In this paper we return to actions of the form (1.1) and reconsider the three-point 
function. The analyses discussed above controlled their calculations by invoking some form 
of slow- roll approximation, typically restricting attention to the lowest power of e = —H/H^ 
(or other quantities of similar magnitude) where H is the Hubble parameter and an overdot 
denotes a time derivative. One may expect this procedure to yield estimates accurate to a 
fractional error of order e, which in some models may be as large as 10~^ to 10~^. Where 
/nl is small, as in canonical scenarios, this is an excellent approximation. However, in 
models where /nl is numerically large, a fractional error of order e may be comparable to 
the resolution of forthcoming data from the Planck survey satellite, especially if the 0(e) 
term enters with a relatively large coefficient. In the equilateral mode, Planck is expected 
to measure the amplitude of three-point correlations with an error bar in the range 25 - 
30, and a future CMB satellite such as CMBPol or CoRE may even achieve A/nl ~ 10 
[53, 54]. Ideally, we would like the theoretical uncertainty in our predictions to fall below 
this threshold. 

Corrections to the power spectrum at subleading order in e are comparatively well- 
understood. For the two-point function of a canonical field, the necessary expressions were 
obtained by Stewart & Lyth [55], who worked to second order in the slow-roll expansion.^ 
Even more precise results, accurate to third order, were given by Gong & Stewart [57, 58]. 
Stewart &; Lyth's calculation was generalized to the noncanonical case by Chen et al.^ [43], 
although third-order results analogous to Gong & Stewart's are not yet available. The second- 
order calculation was systematized by Lidsey et al. [60], who organized their expansion into 
"lowest-order" terms containing the fewest slow-roll parameters, followed by "next-order" 



^Stewart & Lyth were obliged to assume that the slow-roll parameter e — ~H/H^ was small. Grivell & 
Liddle [56] gave a more general construction in which this assumption was not mandatory, but were unable 
to obtain analytic solutions. Their numerical results confirmed that the Stewart-Lyth formulae were valid 
within a small fractional error. 

^Earlier results had been obtained in certain special cases, by Wei et al. [59]. 
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terms containing a single extra parameter, "next-next-order" terms containing two extra 
parameters, and so on. In what follows we adopt this organizational scheme. 

Corrections to /nl have received less attention, but are more likely to be required. The 
prospect of a sizeable "theory error" in a lowest-order slow-roll calculation was identified 
by Chen et al. [43], who computed next-order corrections for P{X,(j)) models. But their 
final expressions were given as quadratures which were not easy to evaluate in closed form. 
How large a correction should we expect? The third-order action derived in Ref. [42] is 
perturbative in the amplitude of fluctuations, but exact in slow-roll quantities: it does not rely 
on an expansion in slow-roll parameters, although it mixes lowest-order and next-order terms. 
Therefore a subset of 0(e) corrections to the vertices of the theory may be obtained without 
calculation, and can be used to estimate the size of the unknown remainder.^ In the case 
of DBI inflation [33, 61], it can be checked that these terms generate a fractional correction 
lOle/7 ~ 14e in the equilateral limit. For the value e ~ 1/20 suggested by Alishahiha, 
Silverstein & Tong [33] this is of order 70%.^ In the absence of further information one 
should expect the unknown terms to be of comparable magnitude but unknown sign, making 
the uncertainty in a purely lowest-order prediction of order 10%-20% even if e is substantially 
smaller than that required for the model of Ref. [33]. In models where /nl is observable, 
but e is not negligible, the lowest-order prediction is likely to survive unaltered only if an 
accidental cancellation removes most of the 0(e) terms. 

There are other reasons to pursue next-order corrections. In models based on (1.1), 
only two "shapes" of bispectrum can be produced with amplitude enhanced by a low sound 
speed. At subleading order, we may expect more shapes to appear — the question is only 
with what amplitude they are produced. If these shapes are distinctive, and the amplitude 
is sufficiently large, it may be possible to detect them directly. Were such a signal to be 
found, with amplitude appropriately correlated to the amplitude generated in each lowest- 
order shape, it would be a striking indication that an action of the form (1.1) was responsible 
for controlling the inflationary fluctuations. 

The calculation presented in this paper enables a precise evaluation of the slow-roll 
corrections, resolving the large theoretical uncertainties and providing detailed information 
about their shapes. For canonical models where |/nl| ^ 1 we verify that the effect is small. 
For noncanonical models we find that the next-order corrections may be enhanced by a 
numerical prefactor, and can be several tens of percent in interesting cases. Although next- 
order calculations are likely to be sufficient for Planck, this suggests that next-next-order 

*These are the next-order components of the gi in Eqs. (3.10)-(3.11), for which it is not necessary to 
compute corrections from the propagator or account for time-dependence of the interaction vertices. If the 
subset is representative of the next-order terms in sign and magnitude then it gives a rough estimate of the 
possible effect. In fact, we will see in §§6-7 that the signs in this subset are probably systematically biased: 
when computed in standard scenarios the remaining terms are of almost exactly the same magnitude but 
opposite sign. But this could not be predicted in advance. 

^Baumann & McAllister later suggested that the Lyth bound [62] placed a limit on e [63]. Lidsey & Huston 
[64] argued that in combination with the large value of e implied by Alishahiha, Silverstein & Tong [33] (see 
also Ref. [65]) this made the "UV" version of the DBI model microscopically unviable. The UV model has 
other difhculties. Bean et al. [66] noted that backreaction could invalidate the probe brane approximation, 
spoiling inflation. For the purpose of making an estimate we are ignoring these details. 
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results could be required to bring the theoretical uncertainty below the observational reso- 
lution of a fourth-generation satellite such as CMBPol or CoRE. If desired, these could be 
computed using the Green's function formalism of Gong &; Stewart [57, 58], although we 
do not attempt this here. We find that P{X^ (p) models can produce a distinctive — but not 
unique — bispectrum shape at next-order. A similar shape arises in certain Galileon models. 

Outline. — This paper is organized as follows. In §2 we discuss the background model and 
specify our version of the slow-roll approximation. In §2.2 we give a short account of the two- 
point function, which is required to study the squeezed limit of the bispectrum. We review 
the derivation of the third-order action in §3, and give an extremely brief introduction to the 
"in-in" formulation for expectation values in §3.1.1. We give a careful discussion of boundary 
terms in §3.1.2. 

Next-order corrections to the bispectrum are calculated in §3.2. Readers familiar with 
existing calculations of the three-point function may wish to skip directly to this section, and 
refer to §§2-3.1.2 for our notation and definitions. We translate the next-order bispectrum 
to a prediction for /nl in §3.4, and study its shape and scale dependence in §§3.5-3.6. We 
identify a new bispectrum shape orthogonal to the lowest-order possibilities, and give cosines 
for the next-order shapes with standard templates. In §3.6 we obtain a general formula for 
the overall running of /nl with scale. 

In §4 we briefly recapitulate the next-order computation for tensor fluctuations. This 
extends the range of observables to include the tensor-to-scalar ratio, r, and the tensor 
spectral index nj, which are helpful in assessing observational signatures. In §§5-6 we apply 
our results to a selection of concrete models. §5 concentrates on canonical single- field inflation 
with an arbitrary potential. The next-order corrections are small, but we are able to verify 
Maldacena's consistency condition to second-order in slow-roll parameters. In §6 we discuss 
models which can imprint a significant nongaussian signature, focusing on DBI inflation and 
/c-inflation. For DBI inflation one can choose a late-time power-law attractor solution, or a 
generic quasi-de Sitter background which depends on the shape of the potential and warp 
factor. In §6.1.1 we calculate the next-order corrections for both cases, obtaining corrections 
to /nl arising from the shape of the potential and warp factor for the first time. We compare 
our results to known exact solutions. In §6.1.2 we apply our calculation to the power-law 
solution of /c-inflation. Finally, we summarize our findings and conclusions in §7. 

A number of appendices collect extra material. We briefly recount the derivation of 
next-order corrections to the propagator in Appendix A. Appendices B-C give mathematical 
details of certain integrals and discuss the role of analytic continuation in obtaining expres- 
sions valid for arbitrary momentum shapes. In Appendix D we provide explicit formulae for 
the quadratures used in Ref. [43] . 

We choose units in which c = h = 1. We use the reduced Planck mass, Mp = (87rG)~^/^, 
usually setting Mp = 1. Spacetime tensors are labelled by Latin indices {a, 6, . . .}, whereas 
purely spatial tensors are labelled by indices {i, j, . . .}. 
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2 Single-field infiation: an overview 



Consider the theory (1.1), with a homogeneous background solution given by the Robertson- 
Walker metric, 

ds"^ = -dt^ + a^{t)dx\ (2.1) 
The corresponding Priedmann equations are 

3H'^ = {2XP^x - P) and 2H + 3^^ = -P, (2.2) 

where H = a/a is the Hubble parameter, an overdot denotes differentiation with respect 
to cosmological time, t, and a comma denotes a partial derivative. The nontrivial kinetic 
structure of P causes fluctuations of the scalar field (p to propagate with phase velocity — 
which we will refer to as the sound speed, Cg — different to unity [67] , 

In the special case of a canonical field, we have P = X/2 — V{(f)) and the fluctuations in (p 
propagate at the speed of light. 

2.1 Fluctuations 

During inflation the universe rapidly isotropizes [68], making (f) spatially homogeneous to a 
good approximation. At the same time, quantum fluctuations generate small perturbations 
whose statistics it is our intention to calculate. Since (p dominates the energy density of the 
universe by assumption, its fluctuations must be communicated to the metric. 

Our freedom to make coordinate redefinitions allows these metric and field fluctuations 
to be studied in a variety of gauges [69, 70]. Without commitment to any particular gauge, 
a perturbed metric can be written in terms of the Arnowitt-Deser-Misner (ADM) lapse and 
shift functions [71], respectively denoted N and Ni, 

ds^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^dt) . (2.4) 

In the absence of perturbations, Eq. (2.4) reduces to the background metric (2.1). The utility 
of this formulation is that N and Ni do not support propagating modes, which are restricted 
to (j) and the spatial metric hij. There are three scalar and two tensor modes, of which one 
scalar can be gauged away by choice of spatial coordinates and another by choice of time. At 
quadratic order the tensor modes decouple from the remaining scalar. The tensors contribute 
at tree-level for n-point functions with n ^ 4 [29, 48], but only at loop-level^ to the three- 
point function. In what follows we work to tree-level and discard tensor modes. Whether or 
not they are retained, however, one uses constraint equations (discussed below) to express 
the lapse and shift algebraically in terms of the propagating modes. ^ 

®The loop expansion in question is generated by insertion of vertices with extra factors of the fluctuations, 
and so is effectively an expansion in powers of ~ 10"^". Conditions under which the "classical" tree graphs 
dominate loops were discussed by Weinberg [72-74] and van der Meulen & Smit [75]. 

'^If desired, the constraints need not be solved explicitly but can be enforced via the use of auxiliary fields, 
yielding an "off-shell" formulation [76]. 
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Where only a single scalar field is present, it is convenient to work on spatial slices of 
uniform field value. In this gauge the propagating scalar mode is carried by hij, and expresses 
local modulations in the expansion a{t), 



h 



(2.5) 



Writing In a{t) = J dN, where dN = H dt and N{t) - N{to) expresses the number of e-folds 
of expansion between times t and to) it follows that C = SN.^ 

Slow-variation parameters. — We work perturbatively in During inflation the universe is 
typically well-approximated by a de Sitter epoch, in which the Hubble parameter is con- 
stant. To express the action (1.1) in terms of C, it is convenient to introduce dimensionless 
parameters which measure the variation of background quantities per Hubble time. We define 



dlnH 
dN 



H 



V 



dine 
dN 



£ 

We' 



and 



din Cg 

dN 



He. 



(2.6) 



Note that e ^ 0, provided H decreases with time. There is no requirement of principle that 
any of e, r/ or s are small, although whenever background quantities are slowly varying we 
may expect e, |s| <C 1. 

Aetion and eonstraints. — In ADM variables, the action (1.1) can be written 



S 



d'^x VhNlR'^^^ + 2P{X,i 



+ 



1 



d^x VhN' 



E' 



(2.7) 



where Eij satisfies 



E: 



{21 



where | denotes a covariant derivative compatible with hij, and indices enclosed in brackets 
(• • • ) are to be symmetrized. The extrinsic curvature of spatial slices, Kij, can be written 
Kij = N-^Eij. 

Varying the action with respect to N and Ni yields constraint equations [42], 



i^(3) ^2P -APxiX + h'^dicpdi 



1 

iV2 



E' 







and 



N 



{Eij — Eh' 



2P 



X 



N 



(2.9a) 



(2.9b) 



Eqs. (2.9a)-(2.9b) can be solved order-by-order. We write N = 1 + a. Likewise, the shift 
vector can be decomposed into its irrotational and solenoidal components, N = diO + j3i, 
where di(3i = 0. We expand a, /3 and 6 perturbatively in powers of C,, writing the terms of 



*Our notation is chosen to coincide with the recent hterature, including Refs. [26, 37, 43, 77-79]. However, 
the uniform-(/) sHcing corresponds to comoving gauge — in which the perturbation in Eq. (2.5) has traditionally 
been denoted TZ [3, 80, 81]. The symbol C, has often been reserved for the corresponding quantity in uniform 
density gauge. The comoving and uniform density slicings agree for adiabatic, superhorizon-scale perturbations 
in single-field inflation [3] , and therefore = 7?. up to an irrelevant sign convention. 
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j^th Qj.(jgj- g^g Q,^^ g^j^d study the three-point correlations it is only necessary to solve 

the constraints to first order [26, 43]. One finds 

ai = -|, /3ii = 0, and = + (2-10) 

The quantity S measures the second X-derivative of P. In what follows it will be helpful to 
define an analogous quantity for the third derivative, denoted A [42], 

E = XP^x + 2X^P,xx = (2.11a) 
A = X'^P^xx + \x^P,xxx. (2.11b) 
2.2 Two-point correlations 

To obtain the two-point statistics of C, one must compute the action (2.7) to second order, 
which was first accomplished by Garriga & Mukhanov [40] . Defining a conformal time variable 
by r = dt/a{t), one finds 

S2 = j d^xdra^zliC'f-clidCf], (2.12) 

where z = e/c^, a prime denotes differentiation with respect to r, and d represents a spatial 
derivative. Although our primary interest lies with P{X, i;A)-type models of the form (1.1), we 
leave intermediate expressions in terms of z and its corresponding slow-variation parameter, 

dlnz i , , 

The last equality applies in a P{X, (j)) model. When calculating three-point functions in §3, 
we will do so for an arbitrary z. In order that the fiuctuations are healthy and not ghostlike 
we must choose z positive. 

To simplify certain intermediate expressions, it will be necessary to have an expression 
for the variation 6S2/SC, 

X = eC-Hx-^d-^^-^, (2.14) 
where we have introduced a variable which satisfies 

d\ = C (2.15) 
The equation of motion follows by setting 882/ 6C, = in (2.14). 

Slow-variation approximation. — Although Eqs. (2.12) and (2.14)-(2.15) are exact at linear 
order in C, it is not known how to solve the equation of motion (2.14) for arbitrary back- 
grounds. Lidsey et al. [60] noted that the time derivative of each slow- variation parameter is 
proportional to a sum of products of slow-variation parameters, and therefore if we assume 

0<e<l, |r/| <Cl, |s| < 1 and |f | < 1, (2.16) 
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and work to first order in these quantities, we may formally treat them as constants. 

Expanding a general background quantity such as H{t) around a reference time t^, one 

finds 

H{t) « H{U) [1 + e,AN,{t) + •••], (2.17) 

where AN^{t) = N{t) — N(t^). Similar formulae apply for the slow-variation parameters e, ry, 
s and V themselves, and their derivatives. Physical quantities do not depend on the arbitrary 
scale t^, which plays a similar role to the arbitrary renormalization scale in quantum field 
theory. Eq. (2.17) yields an approximation to the full time evolution whenever \e-),AN^,{t)\ <^ 
1 [57, 58, 82], but fails no later than ~ l/Si, e-folds after the reference time t^. In order that 
solutions of (2.14) are sufficiently accurate for their intended use — to compute correlation 
functions for t ~ t^, — we must demand that 0(e) quantities are sufficiently small that (2.17) 
applies for at least a few e-folds around the time of horizon crossing, but it is not usually 
necessary to impose more stringent restrictions. This approach will fail if any slow- variation 
quantity becomes temporarily large around the time of horizon exit, which may happen in 
"feature" models [83-87]. 

Working in an arbitrary model, results valid many e-folds after horizon exit typically 
require an improved formulation of perturbation theory obtained by resumming powers of 
AA^ [88-90], for which various formalisms are in use [13, 91-93]. However, as is well-known 
(and we shall see below), this difficulty does not arise for single-field inflation. 

Two-point function. — The time-ordered two-point correlation function is the Feynman propa- 
gator, (T C(r, xi)^(r', X2)) = G{t, t'; |xi — X2I), which depends on the 3-dimensional invariant 
|xi — X2I. In Fourier space G = J d^q {2'iT)~^Gq{T,T')e^'^'^^'^~^^\ which implies 

(T C(ki, r)C(k2, r')) = (27r)35(ki + k2)Gk{r, r'). (2.18) 

The (5-function enforces conservation of three-momentum, and A: = |ki| = |k2|. One finds 

G.(r.r') = ("^M"«^;' '1^7'. (2.19) 

The mode function Uk is a positive frequency solution of (2.14) with 6S2/SC = 0. Invok- 
ing (2.16) and working to first-order in each slow- variation parameter, we find 



(2) 

in which is the Hankel function of the second kind of order u, and zu = e + v/2 + 3,s/2. 
At sufficiently early times, for which |/i:Csr| ^ 1, an oscillator of comoving wavenumber k 
cannot explore the curvature of spacetime and feels itself to be in Minkowski space. In this 
limit (2.20) approaches the corresponding Minkowski wavefunction [94]. 

Power spectrum. — When evaluated at equal times, the two-point function defines a power 
spectrum P{k, r) by the rule 

P{k,T) = Gk{T,T). (2.21) 
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In principle, the power spectrum depends on time. It is conventionally denoted P{k) and 
should not be confused with the Lagrangian P(X,(p). Using (2.19) for r' — )• r, working in 
the limit |A:Csr| — )• and expanding uniformly around a reference time r*, one finds 



^^^^ " AzScl 



2k 

l + 2<{ro^(2-7E-ln— ) 



(2.22) 



The Euler-Mascheroni constant is 7e ~ 0.577. We have introduced a quantity /c^ satisfying 
\ki,Cs-t,Ti,\ = 1, where we recall that is arbitrary and need not be determined by k. Since 
\ki,Csi,Ti,\ « \k^Csi,/ai,Hi,\, we describe as the horizon-crossing time associated with the 
wavenumber k^. Inclusion of next-order effects marginally shifts the time of horizon exit, 
leading to a small mismatch between |A;^Cs^,Tj,| and \ki,Csi,/a^Hi,\.^ In the terminology of 
Lidsey et al. [60], the "lowest-order" result is the coefficient of the square bracket [•••], and 
the "next-order" correction arises from the term it contains of first-order in slow-variation 
parameters. Therefore the lowest-order result can be recovered by setting the square bracket 
to unity. This convention was introduced in Ref. [60] , and when writing explicit expressions 
we adopt it in the remainder of this paper. 

Although expanded around some reference time r^, Eq. (2.22) does not depend on 
AN^, and therefore becomes time-independent once the scale of wavenumber k has passed 
outside the horizon. This is a special property of single-field infiation. Working in classical 
perturbation theory it is known that C becomes constant in the superhorizon limit provided 
the fluctuations are adiabatic and the background solution is an attractor [95-97].^° We are 
not aware of a corresponding theorem for the correlation functions of C, computed according 
to the rules of quantum field theory. However, it appears that in all examples compatible 
with the classical conservation laws discussed in Refs. [95-97], a time-independent limit is 
reached. We will return to this issue in §3.3 below. 

Scale dependence. — The logarithmic term in P{k) indicates that the power spectrum varies 
weakly with scale k, making (2.22) quantitatively reliable only if the reference time k^, is 
chosen sufficiently close to k that | ln(2fc/A:*)| < 1. Defining a "dimensionless" power spectrum 
V by the rule V = k^P{k) /2Tr'^ , the variation of V{k) with scale is conventionally described 
in terms of a spectral index, 

dln-p 

n, - 1 = — — = -2ro^, 2.23 
d m A; 

which is valid to lowest-order provided k-^ w k. Eq. (2.23) is a renormalization group equation 
describing the fiow of V with k, where /3-p = (n^ — 1)V plays the role of the /3-function. An 
expression for Ug valid to next-order can be obtained by setting k = k^, making 'i*r' the 
time of horizon exit of wavenumber k. Having made this choice, the /c-dependence in (2.22) 
appears only through the time of evaluation and is accurate to next-order [55]. We define 

^Chen et al. [43] adopted a definition in terms of |fc*/a*i/*|, making some intermediate expressions different 
in appearance but identical in content. 

^"For canonical inflation this conclusion can be reached without use of the Einstein equations. Recently, 
Naruko & Sasaki argued that the Einstein equations may be required for some types of noncanonical models 
[97]. 
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model lowest-order next-order 



RrbitrRry 


-2e-t 


— 1)^ — 3.Sj, 


— 2e^ -f e^r7^f2 — 


^ /E 


o, 2fc. 
2 In — 


2k 

) + s^t^ (A — 3^t7 — 3 In — 
2A: 

v*u)*(2 - 7e - In — ) 












-f + 


canonical 




- V* 


-2e2+e^7?^(l- 


27E- 


^, 2fc. 
21n — 


2fc 

) +77.e*(2-7E-ln— ) 


P(X,0) 






-2e2+e,?7.(l- 


27E- 
, 2fc. 


2fc. 
21n — 


2fc 

) - s*<*(7e + In — ) 
ft* 












' 3e^s^ - s^ry^ 



Table 1. ?is — 1 at lowest-order and next-order. The first row applies for arbitrary positive, smooth 
z, as explained below Eq. (2.13). 



additional slow-variation parameters, 

and quote results in Table 1 for — 1 at lowest-order and next-order. 



3 Three-point correlations 
3.1 Third-order action 

Three-point statistics can be obtained from the third-order action. This calculation was first 
given in Ref. [42], where the three-point function was obtained under certain hypotheses. 
Chen et al. later computed the full three-point function [43]. After integration by parts, 
using both the background equations of motion (2.2) and the solutions of the constraints 
given in (2.10), we find 

+ d3xdta3|-2-^C(5C)' + 6-^CC'-2^^C'-^52^i9,MiC (3.1) 

where Jq denotes an integral over a formal boundary, whose role we will discuss in more detail 
below. We have temporarily reverted to cosmic time t, rather than the conformal time r. 
The parameters S and A were defined in Eqs. (2.11a) and (2.11b). After further integration 
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by parts, and combining (2.14) and (2.10) with (2.15), one finds^^ 

S, 3 1 /^d». - iSHV + ^(l - |)C(8C)^ - ^«^C(aC)^ - IjCC^ 

+ \ f (l^xdta^\^{e{l-cl)+r]e + e^ + £r]-2es}Q{dCf (3.2) 
+ \{Qe{cI - 1) + 2e2 _ 2er/}CC' 

where / is defined by 

' (3.3) 

Boundary terms. — The boundary terms in Eqs. (3.1)-(3.3) arise from integration by parts 
with respect to time, and were not quoted for the original calculations reported in Refs. [26, 
42, 43]. Adopting a procedure initially used by Maldacena, these calculations discarded all 
boundary terms, retaining only contributions proportional to 882/ liQ in the bulk component 
of (3.2). The 6S2/SC, terms were subtracted by making a field redefinition. 

This procedure can be misleading. The terms proportional to 882/ give no contribu- 
tion to any Feynman graph at any order in perturbation theory, because 6S2/SC is zero by 
construction when evaluated on a propagator. Therefore these terms give nothing whether 
they are subtracted or not. On the other hand, a field redefinition may certainly shift the 
three-point correlation function. Therefore, in general, the subtraction procedure will yield 
correct answers only if this nonzero shift reproduces the contribution of the boundary com- 
ponent in (3.2), which need not be related to /. This argument was first given in Ref. [27], 
and later in more detail in Refs. [90, 98], but was applied to the third-order action for field 



^^In Refs. [42, 43], a further transformation was made to rewrite the terms proportional to rj. Using the 
field equation (2.14) and integrating by parts, these can be consolidated into the coefficient of a new operator 
("^^ — which does not appear in (3.2) — together with corresponding new contributions to / and the boundary 
term. In this paper, we will leave the action as in (3.2) for the following reasons. First, when computing the 
three-point correlation function, the contribution of each operator must be obtained separately. Therefore 
nothing is gained by introducing an extra operator, whose contribution we can avoid calculating by working 
with (3.2). Second, after making the transformation, the contribution from the boundary term is nonzero 
and must be accommodated by making a field redefinition. This redefinition must eventually be reversed to 
obtain the correlation functions of the physical field (. If we leave the action as in (3.2) then it transpires 
that no field redefinition is necessary. 
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fluctuations in the spatially flat gauge. In this gauge only a few integrations by parts are 
required. The boundary term is not complicated and the subtraction procedure works as 
intended. In the present case, however, it appears impossible that the subtraction procedure 
could be correct, because the boundary term contains operators such as which are not 
present in /. Indeed, because C approaches a constant at late times, the term appar- 
ently leads to a catastrophic divergence which should manifest itself as a rapidly evolving 
contribution to the three-point function outside the horizon. 

This potential problem can be seen most clearly after making the redefinition Q ^ — f 
under which the quadratic action transforms according to 

S2[Q\ ^ -2 I a^^nf - [ d^xdr (3.4) 

The bulk term proportional to 6S2/5C, disappears by construction. After the transformation, 
the boundary term becomes 



53 2 \ I A - isH,^ + ^(1 - |).(a.)^ - ^avo.) 



(3.5) 



in which x is to be interpreted as a function of vr [cf. Eq. (2.15)]. 

Eq. (3.5) is not zero. To satisfy ourselves that it does not spoil the conclusions of 
Refs. [26, 42, 43], we must determine how it contributes to the three-point correlation func- 
tion. Before doing so, we briefly describe the in-in formalism which is required. Readers 
familiar with this technique may wish to skip to §3.1.2. 

3.1.1 Schwinger's in in formalism 

The correlation functions of interest are equal time expectation values taken in the state 
corresponding to the vacuum at past infinity. At later times, persistent nontrivial correlations 
exist owing to gravitational effects associated with the time-dependent background of de 
Sitter. 

Feynman's path integral computes the overlap between two states separated by a finite 
time interval, which is taken to infinity in scattering calculations. Schwinger obtained ex- 
pectation values at a finite time t^, by inserting a complete set of states \i,tf) at an arbitrary 
time tf t^:, 

(in|0(t,)|in) = '^{m\i,tf){i,tf\OiU)\m). (3.6) 

i 

where |in) is the "in" vacuum in which one wishes to compute the expectation value and O 
is an arbitrary local functional. Our choice oi tf is irrelevant. Choosing a basis of energy 
eigenstates, \i,tf) = e~^^*^^f~^*^\i,t^,) where Ei is the energy of the state \i). This phase 
cancels in (3.6).^^ Expressing each overlap as a Feynman path integral, we obtain [99-102] 

(in|0(t,)|in) = / [d0+ d,/._] 0{U) e'^['^+l-^^['^-l5[0+(t,) - </>-(*.)]. (3.7) 



^■^One could just as well insert a complete set of states at an arbitrary time t < t», but the resulting overlap 
would not be expressible in terms of a path integral. Had we retained tf > t,, the resulting contributions 
would have yielded only cancelling phases in (3.7). 
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The (5- function restricts the domain of integration to fields cj)j^ and (/>_ which agree at time t^, , 
but are unrestricted at past infinity. The requisite overlap of an arbitrary field configuration 
with the vacuum is obtained by deforming the contour of time integration. Cosmological 
applications of Schwinger's formulation were considered by Jordan [103] and by Calzetta 
& Hu [104], to which we refer for further details. Applications to inflationary correlation 
functions were discussed by Weinberg [73, 74] and have been reviewed elsewhere [105, 106]. 

3.1.2 Removal of boundary terms 

If boundary terms are present they appear in (3.7) as part of the action S with support at 
past infinity and t = . The deformed contour of integration kills any contribution from past 
infinity, leaving a boundary term evaluated precisely at t*, where the 5- function constrains 
the fields to agree. Therefore, at least for O containing fields but not time derivatives of 
fields, any boundary operators which do not involve time derivatives produce only a phase 
which cancels between the + and — contours. This cancellation is a special property of 
the in-in formulation: it would not occur when calculating in-out amplitudes, for which 
the uncancelled boundary term would diverge near future infinity. Rapid oscillations of e'*^ 
induced by this divergence would damp the path integral, yielding an amplitude for any 
scattering process which is formally zero. This can be understood as a consequence of the 
lack of an S-matrix in de Sitter space [107]. 

In virtue of this cancellation we may disregard the first three operators in the boundary 
part of (3.5). However, the (5- function in (3.7) in no way requires that time derivatives of 
the + and — fields are related at t = t^. Therefore operators involving time derivatives need 
not reduce to cancelling phases. To understand their significance we subtract them using a 
further field redefinition. 

Inspection of Eq. (3.5) shows that the time-derivative terms are of the schematic form 
TTvr^, and therefore lead to a field redefinition of the form ^ — >• vr + tttt. We now argue that 
boundary operators with two or more time derivatives are irrelevant on superhorizon scales. 
Using the schematic field redefinition, the three-point correlation functions of C and vr are 
related by {('^) = {ir^) + 3{'ir'^) {tttt) plus higher-order contributions. However, Eq. (2.22) 
implies (tttt) — )• on superhorizon scales, and therefore (C^) = (tt^) up to a decaying mode. 
This field redefinition will inevitably produce bulk terms proportional to 6S2/SC, but we have 
already seen that these do not contribute to Feynman diagrams at any order in perturba- 
tion theory. Therefore, on superhorizon scales, the correlation functions of the original and 
redefined fields agree. It follows that after subtraction by a field redefinition, the unwanted 
boundary terms in (3.5) are irrelevant and can be ignored. Equivalently, one may confirm 
this conclusion by checking that operators with two or more time derivatives give conver- 
gent contributions to the boundary action at late times. Similar arguments apply for any 
higher-derivative combination. 

We conclude that the only non-negligible field redefinitions are of the schematic form 
^ — )■ tt+tt'^, which arise from boundary operators containing a single time derivative. Eq. (3.5) 
contains no such operators. However, had we made a further transformation to consolidate 

^^Spatial derivatives play no role in this argument, which therefore applies to the entire first line of (3.3). 
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the rj dependence, a contribution of this form would have appeared in / and the boundary 
action. It can be checked that this single-derivative term would be correctly subtracted by 
Maldacena's procedure, and in this case the subtraction method applied in Refs. [26, 42, 43] 
yields the correct answer. However, in theories where single-derivative terms already appear 
in the boundary component of (3.1) there seems no guarantee it will continue to do so. 

3.2 The bispectrum beyond lowest-order 

We define the bispectrum, B, in terms of the three-point function. 



Observational constraints are typically quoted in terms of the reduced bispectrum, /nl, which 
satisfies [8, 13] 



Current constraints on /nl in the simplest inflationary models have been discussed by Sen- 
atore et al. [108]. 

Next-order corrections. — We are now in a position to compute the three-point function of (3.2) 
to next-order. The arguments of the previous section show that the boundary action and 
6S2/SC contributions can be discarded. Next-order terms arise in the remaining operators 
from a variety of sources. First, the coefficients of each vertex in (3.2) contain a mixture of 
lowest-order and next-order contributions. Second, the lowest-order part of each vertex is a 
time-dependent quantity which must be expanded around a reference time, as in (2.17), pro- 
ducing next-order terms. Third, there are next-order corrections to the propagator, obtained 
by expanding (2.19)-(2.20) in the neighbourhood of the chosen reference time. Propagator 
corrections appear on both the external and internal legs of the diagram. 

Reference time, factorization scale. — To proceed, we must choose a reference point around 
which to expand time-dependent quantities. Consider an arbitrary correlation function of 
fields C(kj)- Whatever our choice of r^, the result (2.22) for the power spectrum shows that 
we must expect logarithms of the form In ki/k^ which account for the difference in time of 
horizon exit between the mode ki and the reference wavenumber k^,. To obtain a reliable 
answer we should attempt to minimize these logarithms. 

If all fields participating in the correlation function carry momenta of approximately 
common magnitude ki ~ k — described as the "equilateral limit" — the logarithm will be small 
when fc* ~ k. In this case, naive perturbation theory is not spoiled by the appearance of large 
logarithms. In the opposite limit, one or more fields have "soft" momenta of order fciR which 
are much smaller than the remaining "hard" momenta of order kjjy. When ki^i/kuY — >■ it 
will not be possible to find a choice of k^, which keeps all logarithms small and the calculation 
passes outside the validity of ordinary perturbation theory. We have encountered the problem 
of large logarithms which led to the renormalization group of Gell-Mann & Low [109]. 

In the study of inflationary correlation functions, configurations mixing hard and soft 
momenta with /cjr <^ fcuv are referred to as "squeezed," and are of significance because 
they dominate the bispectrum for canonical inflation [26]. In principle one could study the 



(C(ki)C(k2)C(k3)) = {27rfS{ki + k2 + ks)B{ki, k2, ks). 





(3.9) 
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behaviour of a correlation function as its momenta are squeezed by setting up an appropriate 
renormalization group analysis [110]. But this is more complicated than necessary. Malda- 
cena argued that, as the momentum carried by one operator becomes soft, the three-point 
function would factorize: it can be written as a "hard subprocess," described by the two-point 
correlation between the remaining hard operators on a background created by the soft opera- 
tor [26]. Factorization of this kind is typical in the infrared dynamics of gauge theories such as 
QCD, where it plays an important role in extracting observational predictions. The various 
factorization theorems for QCD correlation functions have been comprehensively reviewed by 
Collins, Soper &: Sterman [111].^^ Maldacena's argument was later generalized by Creminelli 
et al. [112]. The factorization property can be exhibited by an explicit decomposition of the 
field into hard and soft modes [113, 114]. 

Because the squeezed limit can be described by Maldacena's method, the outcome of this 
discussion is that the reference scale should usually be chosen to minimize the logarithms 
when all momenta are comparable. In the remainder of this paper we quote results for 
arbitrary fc*, but frequently adopt the symmetric choice k^, = ki + k2 + ks where numerical 
results are required. Having done so, we will be formally unable to describe the squeezed 
limit. Nevertheless, because there is no other scale in the problem, our results must be 
compatible with the onset of factorization in appropriate circumstances — a property usually 
referred to as Maldacena's consistency relation. We will see below that this constitutes 
a nontrivial check on the correctness of our calculation; see also Renaux-Petel [114] for a 
recent discussion of Maldacena's condition in the case of P{X, (p) models. 

Operators. — To simplify our notation, we rewrite the cubic action (3.2) as 



but our calculation will apply for arbitrary gi. Although (" is dimensionless, it is helpful for 
power-counting purposes to think of it as a field of engineering dimension [mass], obtained 
after division by the Hubble rate H. In this counting scheme, the ^'^ operator is dimension-6, 

^■^The background created by soft modes is typically described by some version of the DGLAP (or Altarelli- 
Parisi) equation. A similar phenomenon seems to occur in the inflationary case [IfO]. Equally, the sep- 
arate universe method can be thought of as a factorization theorem for secular time-dependent logarithms 
~ ln|fcCi,r|. The 5N rules which translate correlation functions of the field perturbations into correlation func- 
tions of are an important special case. In this sense, factorization is as important in extracting observable 
quantities for inflation as it is for QCD. 

^^In the gauge theory language discussed above, the scale fc* can be thought of as the factorization scale. 
Operators carrying momentum fc <C fc* should not be included as part of the hard subprocess, but factorized 
into the background. 



^3 




(3.10) 



In a P{X, (j)) model the interaction vertices are 




(3.11) 
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whereas the remaining four operators are dimension-5. At low energies one would naively 
expect the dimension-6 operator to be irrelevant in comparison to those of dimension-5. 
However, the dimension-5 operators are suppressed by the scale H making all contributions 
equally relevant. This manifests itself as an extra power of H in the denominator of gi. 

The vertex factors gi are themselves time-dependent background quantities. We define 
slow-variation parameters hi which measure their rate of change per e-fold, 



and take these to be 0(e) in the slow- variation approximation. 
3.3 Three-point correlations 

We use these conventions to compute the next-order bispectra for each operator in (3.10). 
The resulting three-point functions are complicated objects, and when quoting their values it 
is helpful to adopt an organizing principle. We divide the possible contributions into broadly 
similar classes. In the first class, labelled 'a,' we collect (i) the lowest-order bispectrum; 

(ii) effects arising from corrections to the wavefunctions associated with external lines; and 

(iii) effects arising from the vertex corrections. In the second class, labelled '6,' we restrict 
attention to effects arising from wavefunctions associated with internal lines. These are 
qualitatively different in character because wavefunctions associated with the internal lines 
are integrated over time. Adapting terminology from particle physics, we occasionally refer 
to the lowest-order bispectrum as the "LO" part, and the next-order piece as the "NLO" 
part. 

Large logarithms, infrared singularities. — The computation of inflationary n-point functions 
has been reviewed by Chen [115] and Koyama [106]. At least three species of large logarithms 
appear, disrupting ordinary perturbation theory. We carefully track the contribution from 
each species. The most familiar types — already encountered in the two-point function — 
measure time- and scale-dependence. A third type of large logarithm is associated with the 
far infrared limit A;ii^/A;uv ~^ discussed in §3.2. This is Maldacena's "squeezed" limit, 
discussed in §3.2, in which the behaviour of the three-point function obeys a factorization 
principle. We will show that the various large logarithms arrange themselves in such a way 
that they can be absorbed into the scale-dependence of background quantities. 

Time-dependent logarithms appear after expanding background quantities near a fixed 
reference scale, as in (2.17), where at conformal time r we have A'^^ = ln|A;^Csr|. In §2.2 
we explained that the correlation functions of C are expected to become time-independent 
outside the horizon. Therefore one should expect all In r dependence to disappear. Some A''*- 
type logarithms cancel among themselves but others cancel with time-dependent logarithms 
arising from wavefunction corrections associated with internal lines. The internal lines are 
aware only of the intrinsic geometrical scale kt = ki + k2 + k^ and cannot depend on the 
arbitrary reference scale k^, so the outcome of such a cancellation leaves a residue of the 
form Inkt/ki,. These are scaling logarithms, entirely analogous to the logarithm of (2.22), 
describing variation of the three-point function with the geometrical scale kf. Scale logarithms 
can also occur in the form In/cj/A;^. 



hi = 



Hg^ 
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The third species of logarithm takes the form Inki/kt- Each side of the triangle must 
scale linearly with the perimeter, so despite appearances these have no dependence on kt — 
they are unaffected by rigid rescalings of the momentum triangle (cf. Eqs. (3.23a)-(3.23c) 
below). We describe them as 'purely' shape dependent. The 'pure' shape logarithms become 
large in the squeezed limit ki/kt — )■ 0. 

a-type bispectrum. — Collecting the o-type contributions to the bispectrum, we find 
_ Ht T°(fci) f k,k2ks 



+ 2F"(A:i)e* In ^ + Wiki^hi^ In ^ 

^ ^ A:* ^ k^ (3.13) 



+ X''{ki){l + 3^^) + 2y"(fci)e^ + Z''{ki)hi 
+ cyclic permutations. 

The coefficients r'^(fei), C/"(fei), l/"(A:i), ^{ki), X"(A;i), ^^(/ci) and Z"(A:i) are functions 
of all three momenta kj and are symmetric under the exchange ^2 -H- k^. We adopt the 
convention, used through the remainder of this paper, of writing only the asymmetric mo- 
mentum explicitly. The notation 'cyclic permutations' denotes addition of the preceding term 
under cyclic permutations of the kj. The result is symmetric under interchange of any two 
momenta. 

We give explicit expressions for the coefficient functions in Table 2. The quantity ii^ is a 
combination of slow- variation parameters, E = w[2 — 7e — ln2) — e — s, and also appears in 
the power spectrum (2.22). The term proportional to X"'{ki) includes the entire lowest-order 
bispectrum. 

b-type bispectrum. — The 6-type bispectrum must be added to the a-type terms. It has no 
lowest-order contributions, and can be written 



B 



+ (^X''J2{ki) + Y^'kf In + + cs,k^ Re (J.) ^ ^^'^^^ 
+ cyclic permutations. 

The same convention applies to the arguments of the coefficient functions T^{ki), U^{ki), 
Vf'{ki), W''{ki), X^ Y'' and Z^. We give exphcit expressions in Table 3. 

Eq. (3.14) depends on three logarithmic functions Jj (which are not Bessel functions) 
defined by 

2k- 

AJo{ki) = ln-^, (3.15a) 
kt 

2k- 

T?fJi(A;i) = i?, + ln-^, (3.15b) 
kt 

2k- 

dp2{ki) = ^i{2 + di) + 2\u-^, (3.15c) 

h 
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y(k2 •k3) 


fc?(k2-k3) 


Orl* ; 


X (fci + kt) 
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1 






3kt — ki 


1 




1 


1 




^(^^2 ^2 ^lfc2fc3^ 


3kt — ki 


1 




1 


1 


3cL 


(a^2 ^2^fclfc2fc3^ 


3kt — ki 


1 




1 


1 


cL( 


^2^fclfc2fc3^ 


3kt — ki 


1 




1 

7E-2 


^ ^* 








(3fct - fci)7E 


^ ^* 


fci + fct 


f,2 7^2 kik2k3\ 




+ 2kt 


fci + kt 



Z^ik,) 7E 



7e 



ki 



3ci 



fci + kt 



+ (1-7e)U: 



2 fcifc2fc3 



(3fct - fci)7E ki 
7e 

+ fci - fct 



fci + kt 



Table 2. Coefficients of tfie leading order bispectrum. = fcifc2 + A;ifc3 + k2kj,. 

where i?j = 1 — 2ki/kt. These exhaust the 'pure' shape logarithms of the form Inki/kt, 
discussed above, which appear only in the Jj. There is an obvious logarithmic divergence 
in the squeezed limit ki — >• 0, which we will show to be responsible for factorization of the 
correlation function. There is potentially a power-law divergence in the limit kt — t- 2ki. This 
is also a squeezed limit — in which the i*^ side stays fixed while a different momentum goes 
to zero. In this limit T?j — >• 0, making the Jj naively divergent. If present, such power-law 
divergences would be puzzling. However, it can be checked that — in combination with the 
logarithm — each Jj is finite. This infrared-safe behaviour relies on a resummation procedure 
which is discussed in more detail in Appendix B. 
The function satisfies 



-ii + 5i 272 + 3^2 673 -m (^3 2474 + 50(^4 

fc2 + 1,3 + l4 

i t t 

, 1 ^* 1 /^A '^l 0*^2 , „(53 , ^. ^4 



(3.16) 



This function is discussed in Appendix C. The coefficients 70, 71, 72, 73, 60, Si, 62 and S3 
depend on the operator under consideration. We quote values for each operator in Table 4. 
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operator 
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— i-1 
ft,l 






ft,l 




^ (k-, — k^) 
2\'^3 K2) 


— i-1 








kt — 2k 








1 
















2 












kf[2TUu - 3 Re(^o^) + 37_Bn7i^] 







Table 3. CoefRcients of the subleading corrections to the bispectrum. ^ kik2 + kik^ + fc2fc3- 

The operators CC'^j C'9jCd~'^djC, and d'^C^{d~'^djC,')'^ are all dimension-5, and differ only 
in the arrangement of spatial gradients. For arbitrary shapes their three-point functions will 
not coincide, but for equilateral triangles the arrangement of gradients is irrelevant and the 
resulting /nl should agree. This will represent a minimal check of our expressions. We will 
carry out further checks in §3.6 and §§5-6. 

A subset of these terms were calculated by Chen et al. [43]. Our calculations exhibit 
two principal differences. First, Chen et al. worked to fixed order in slow-roll quantities, 
keeping terms of 0(e) only. In a model where <C 1 this gives the next-order corrections. 
However, in a model where ~ 1 the leading terms are themselves 0(e) and the formulae of 
Chen et al. reduce to these leading contributions. In our calculation, we work uniformly to 
next-order rather than a fixed order in powers of e. When ^ 1 our next-order corrections 
are 0(e), and we have verified that they agree with those computed in Ref. [43]. (We give 
more details of the relation between our calculations in Appendix D.) When ~ 1 the 
next-order corrections are O(e^). These were not included in the formulae of Ref. [43]. 

Second, we retain a fioating reference scale /c^. In Appendix B^^ of Ref. [43] this was cho- 

^® Slightly different conventions for fc* were used elsewhere in Ref. [43]. 
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operator 
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/io* + 2s* - 2^1* 


7,2 1 ly, , , ( 1^ 1 ^ 




/io* + 2s* - 2/ii* 
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-Si,k\{ki + fca) 

- ptl*fclfc2fc3 


/ctS* - 3fc2S* 
+ 3/ii*fc2 


fcts* - 3fcis* 
+ 3/ii*fci 


72 






-Si,k\kikj, 


k2ktSi, 


fclfctS* 




73 


I, 

r 










-50 




37Z71* - 4s* 




3roi* - 4s* 


3n7i* - 4s* 


6i 






Si,kl{k2 + ks) 
+ mii,kik2k3 


-fctS* + 5fc2S* 
- 3^2^!* 


-fcts* + 5fcis* 
- 3fcin7i* 




tui* — 2s* 

"-Sit 




Si,klk2k3 


-k2ktSi, 


-fclfctS* 















Table 4. Coefficients appearing in the function J' for each operators. Note that the 7i contain complex 
numbers. The imaginary part is cancelled on addition of the + and — Feynman diagrams, and only 
the real part of these coefficients contribute. In an intermediate step for the three-point function 
of C(9C)^, the cancellation of power-law divergences in the conformal time r (which is required by 
Weinberg's theorem [73]) depends on a real contribution generated from the product of two imaginary 
terms. 



sen to be fc* = fcf. Retaining this scale allows us to extract the scale- and shape-dependence 
of /nl (§§3.5-3.6). 

3.4 Formulae for /nl 

The individual bispectra, with their detailed shape-dependence, are the principal observable 
objects. However, for simple model comparisons it is helpful to have an explicit expression 
for the nonlinearity parameter /nl defined in Eq. (3.9). Accounting for scale-dependent 
logarithms present in the power spectrum, one finds 

^5/4z*c|*y B{h,k2,k^)U^kf 

Q \ E.kf{l + AE^-2^uln{k7^k;'U,k,}y ^ 

This expression is to be expanded uniformly to 0(e) in slow- variation parameters. 
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There is another reason to study /nl- We have explained that large logarithms of 
the form In ki/ki, or In ki/kt are to be expected in the squeezed limit /cj — )• 0, describing 
variation of the bispectrum with shape. The power spectrum P{k) contains similar large 
logarithms. Since copies of the power spectrum must be factored out to obtain /nl, one may 
expect it to be more regular in the squeezed limit. Indeed, a stronger statement is possible. 
Partitioning the momenta into a single soft mode of order km and two hard modes of order 
feuV) Maldacena's consistency condition requires [26] 



/NL^-^(n.-l)|fcu^, (3.18) 



as feiR 0, where the right-hand side is to be evaluated at horizon exit for the mode of 
wavenumber fcuv- Eq. (3.18) is finite and independent of any logarithms associated with 
the limit km — )■ 0, which is why this behaviour is described as factorization. It imposes the 
nontrivial requirement that all large logarithms can be absorbed into P(fciR). Such logarithms 
are subtracted by the denominator of (3.17), making /nl finite. 

For each operator i, we write the corresponding /nl as /nl^ and quote it in the form 

hu = /nlIjO [l + l^h\ihi* + K^\iV^ + KsliS* + l^e\i^*] ■ (3.19) 

In Tables 5 and 6 we give explicit expressions for the coefficient functions /nlIjo 
in the case of equilateral^^ and squeezed triangles. Table 6 confirms that (3.19) is finite 
in the squeezed limit, as required. In the equilateral case, we find that the operators CC'^j 
C'djCd~'^djC and d'^C{d~'^djC')'^ agree, for the reasons explained above. 

3.5 Shape dependence 

Stewart &: Lyth's interest in next-order corrections to the power spectrum lay in an accurate 
estimate of its amplitude. In comparison, next-order corrections to the bispectrum could 
be relevant in at least two ways. First, they could change the amplitude of three-point 
correlations, as for the power spectrum. Second, they could lead to the appearance of new 
"shapes," by which is meant the momentum dependence of B{ki, k2, k^) [118], defined in (3.8). 
In principle, both these effects are measurable. 

3.5.1 Inner product and cosine 

Babich et al. introduced a formal "cosine" which may be used as a measure of similarity in 
shape between different bispectra [118]. Adopting Eq. (2.21) for the power spectrum P{k), 
one defines an inner product between two bispectra Bi and B2 as 



R R _ Bi{ki,k2,k3)B2{ki,k2,k3) 

'''■^-=2^ Pik,)P{k2)Piks^ ' ^'-^'^ 



triangles 



The quoted quantity is fm.{k,k,k), which is not the same object as /J^l"' f^i" which constraints are 
typically derived from data [6, 108]. To obtain one should take an appropriately normalized inner 

product (see §3.5 for a simple example) between the full next-order bispectrum and the equilateral template 
[116, 117]. At this level of precision, it may even be desirable to include experiment-dependent information 
in the inner product. 
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where the sum is to be taken over all triangular configurations of the kj. The cosine between 
Bi and B2 is 

B B 

cos(Bi,B2) = ^ -77:7. (3.21) 

These expressions require some care. In certain cases the result may be infinite, requiring 
the summation to be regulated. 

Inner product. — We define the sum over triangles as an integral over triangular configurations 
in a flat measure, so ^ — )• / d^ki d^k2 d^k^ 5(ki+k2+k3). It is sometimes useful to introduce 
a more complicated measure, perhaps to model observational effects [119]. In this paper we 
retain the flat measure for simplicity. The (^-function can be integrated out immediately, 
leaving a space parametrized by two vectors forming a planar triangle which we choose to 
be ki and k2. The triangle is invariant under a group S0(2) x U(l), representing arbitrary 
rotations of ki combined with azimuthal rotations of k2; these change our representation of 
the triangle but not its intrinsic geometry. The volume of this group may be factored out of 
the measure and discarded. Reintroducing k^, in favour of the remaining angular integration, 
we conclude 

R R _ f (-rr , .j^ Bi{ki,k2,k3)B2{ki,k2,k3) 

B,-B2-J { ii k. dh) • (3.22) 

The ki can be parametrized geometrically in terms of the perimeter, fcj, and two dimensionless 
ratios. We adopt the parametrization of Fergusson &: Shellard [119], 

fci = ^(l + a + /3) (3.23a) 

k2 = j{l-a + P) (3.23b) 

A;3 = |(l-/3), (3.23c) 

where 0^/3^1 and /3 — 1 ^ a ^ 1 — /3. The measure dki dk2 dfcs is proportional to 
kldktdadfi. Also, on dimensional grounds, each bispectrum Bi scales like Bik^^, where 
Bi is dimensionless, and each power spectrum P scales like Pk^^ where P is dimensionless. 
In the special case of scale-invariance, P is constant and the Bi depend only on a and j3. 
Therefore 

Bi-B2 = N j dad/3(l-/3)(l + a + ^)(l-a + /3)Si(a,/3)52(a,/3), (3.24) 

where is a harmless infinite normalization which can be divided out. With this under- 
standing we use (3.24) to determine the cosine of Eq. (3.21). In practice, our bispectra are 
not scale invariant and therefore (3.24) does not strictly apply. However, the violations of 
scale invariance (to be studied in §3.6 below) are small. 

Divergences. — Eq. (3.24) may be infinite. For example, the well-studied local bispectrum 
diverges like (1 + a + /3)~2 or (1 — a + f3)~'^ in the limit /3 — )• 0, a — )• ±1, or like (1 — 
/3)~2 in the limit /3 — ?■ 1, a — [118]. These correspond to the squeezed limits discussed 
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in §3.2. Eq. (3.24) therefore exhibits power-law divergences on the boundaries of the region 
of integration, and in such cases the integral must be regulated to obtain a finite answer. For 
simplicity, we adopt a sharp cutoff which requires ki/kt > Smm- As Smm — ^ the cosine (3.21) 
may converge to a nonzero limit if Bi ■ B2, Bi ■ Bi and B2 ■ B2 diverge at the same rate. 
Otherwise, except in finely-tuned cases, it converges to zero. 

For this reason, where divergences exist, the value assigned to cos(-Bi,i?2) is largely 
a matter of convention. However, to resolve the practical question of whether two shapes 
can be distinguished by observation it should be remembered that experiments cannot mea- 
sure arbitrarily small wavenumbers. Therefore their ability to distinguish shapes peaking 
in the squeezed limit is limited. In this case, to obtain accurate forecasts of what can be 
distinguished, one should restore the /cf-dependence in (3.24) and restrict the integration to 
observable wavenumbers, yielding a manifestly finite answer [119, 120]. 

3.5.2 Bispectrum shapes from slow-variation parameters 

In a model with arbitrary gi, the bispectrum is a linear combination of the shapes produced 
by the five operators in (3.10). Of these, CC'^ ^-nd C(f^C)^ ^-re predominantly correlated with 
the local template and the remainder correlate strongly with the equilateral template. The 
C'^ operator has some overlap with the enfolded template, yielding a cosine of order 0.75. 
With generic values of the slow-variation parameters the situation at next-order is similar, 
and each next-order shape is largely correlated with its parent lowest-order shape. 

Lowest-order shapes. — A P{X, (p) is not generic in this sense, but imposes strong correlations 
among the Qi. At lowest order and 55 do not contribute. We focus on a model with small 
sound speed, in which next-order corrections are most likely to be observable, and retain 
only contributions enhanced by c~'^. The remaining three operators organize themselves into 
a family of shapes of the form Si + 01S2, where ^2 arises only from C,'^ but Si is a linear 
combination of the shapes produced by C,'^, CC''^ and ({dC,)'^. The parameter a is the enhanced 
part of A/S, that is 



In the DBI model a = 1/2. We plot the shapes 5*1 and 52 in Table 9. Note that although 
Si involves a linear combination of the local-shape operators CC'^ arid C(c^C)^) the P{X,cj)) 
Lagrangian correlates their amplitudes in such a way that there is no divergence in the 
squeezed limit. Both 5i and 52 are strongly correlated with the equilateral template. They 
are similar to the Mi- and M2-shapes studied in a Galileon theory by Creminelli et al. [37]. 

Next-order shapes}^ — At next-order, more shapes are available. Naively, the family of en- 
hanced bispectra is labelled by e, r/, s and also £ (following Chen et al. we define i = X/HX 
[43]). In practice there is some degeneracy, because the shapes corresponding to these in- 
dependent parameters may be strongly correlated. We will see these degeneracies emerge 
naturally from our analysis. 

^*We thank Xingang Chen and Sebastien Renaux-Petel for helpful discussions relating to the material in 
this section. 




(3.25) 
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The Cg ^-enhanced next-order shape can be written as a Unear-combination of shapes 
proportional to the e, rj, s and i parameters, modulated by a, 



eSs + rjSn + sSs + a [eS'^ + rjS'^ + sS'g + iS'^) . 



(3.26) 



We give overlap cosines of the Si , 5- with the standard templates in Table 7 and plot their 
shapes in Table 10. Because the cutoff dependence complicates comparison between different 
analyses we list the cosines between templates in Table 8, computed using the same conven- 
tions. Generally speaking, these shapes have strong overlaps with the equilateral template. 
However, two are quite different in appearance and have a slightly smaller cosine ~ 0.85 
with this mode: these are S'^ and Sg- We fix two coefficients in (3.26) by choosing a linear 
combination orthogonal to both Si and 82- Without loss of generality we can choose these 
to be r] and s. We find the required combination to be approximately 



It is possible this procedure is stronger than necessary. Both 5*1 and 52 are correlated with 
the equilateral template, and it may be sufficient to find a linear combination orthogonal 
to that. In what follows, however, we insist on orthogonality with Si and ^2 and defer 
generalizations to future work. For certain values of a the denominator of both r] and s 
may simultaneously vanish, making the required r] and s very large. This implies that, near 
these values of a, no shape orthogonal to both 5*1 and S2 can be found within the validity 
of next-order perturbation theory. Therefore we restrict attention to those a which allow 
acceptably small rj and s. 

This process leaves two linear combinations proportional to e and £. In principle these 
can be diagonalized, yielding a pair of shapes orthogonal to each other and {Si, S2}. However, 
the 2x2 matrix of inner products between these linear combinations is degenerate. Therefore, 
only one member of this pair is physical and can be realized in a P{X, tp) model. The other 
is not: it has zero inner product with (3.26), and is impossible to realize because of enforced 
correlations between coefficients. We denote the physical orthogonal combination O. It 
has a vanishing component proportional to 1. This was expected, because the shape S*^ is 
the same as 5*2. For this reason, Chen et al. [43] absorbed £ into a redefined A/S. It is 
indistinguishable from the lowest-order prediction and could never be observed separately, 
which is the origin of the degeneracy. We could have arrived at the same O by excluding S'^ 
from (3.26). Demanding the inner product with Si and 5*2 be zero reproduces the physical 
linear combination obtained from diagonalization. 

We plot the shape of O in Table 11. Its dependence on a is modest. As a function of the 
ki there are multiple peaks, and therefore O is not maximized on a unique type of triangle. 
In Table 13 we give the overlap cosine with common templates. The lowest-order shapes Si 
and ^2 are strongly correlated with the equilateral template, and since O is orthogonal to 
these by construction it also has small cosines with the equilateral template, of order 10~^. 
There is a moderate cosine with the local template of order ~ 0.3 - 0.4. The precise value 



s 



0.12a£(a + 0.72)(a + 1.82) - 0.88£(q - 9.15)(« - 0.22)(a + 1.82) 

(a - 10.24)(a - 0.23)(a + 1.82) 
a£(3.88 - 0.12a) - 1.12£(a - 8.51)(a - 0.08) 
(a - 10.24)(a-0.23) 



(3.27b) 



(3.27a) 



- 24 - 



depends on our choice of 6nim, but the dependence is not dramatic. In Table 13 we have 
used our convention 6ram = 10~^. For Smm = 10~^ the local cosines change by roughly 25%. 
Overlaps with the remaining templates are stable under changes of Suim- There is a cosine of 
order 0.35 - 0.40 with the orthogonal template, and of order 0.30 - 0.35 with the enfolded 
template. We conclude that O is not strongly correlated with any of the standard templates 
used in CMB analysis. To find subleading effects in the data, it will probably be necessary 
to develop a dedicated template for the purpose. 

The 0-shape is very similar to a highly orthogonal shape constructed by Creminelli 
et al. [37] in a Galileon model, although O contains marginally more fine structure. For 
comparison, we plot the Creminelli et al. shape in Table 12 and include its cosine with O in 
Table 13. For varying a we find a cosine in the range 0.8 - 0.9, which indicates it would be 
difficult to distinguish these shapes observationally. In particular, even if a bispectrum with 
this shape were to be detected, further information would be required to distinguish between 
candidate P{X, (j)) or Galileon models for its origin. 

3.6 Scale dependence 

In this section we use the logarithms In ki/k^, and In kt/k^, to study the scale-dependence of 
the three-point function. In the squeezed limit this is determined by (3.18). The only scale 
which survives is the common hard momentum /cuv ) and the variation of /nl with this scale 
is determined by the variation of — 1. This is typically called the running of the scalar 
spectral index [122], and leads to a further consistency relation inherited from Maldacena's — 
and, in general, a hierarchy of such consistency equations generated by taking an arbitrary 
number of derivatives. In the case of single-field canonical inflation, discussed in §5 below, 
we are able to verify this explicitly. 

Away from the squeezed limit, deformations of the momentum triangle may change 
either its shape or scale. Scale dependence occurs in even the simplest models for the same 
reason that the spectrum V and spectral index Ug depend on scale [13, 26, 27, 42, 43]. Chen 
introduced a 'tilt,' Uf^^, defined by^^ [123] 

For a fixed triangular shape, this measures changes in /nl as the perimeter varies. Scale 
dependence of this type was subsequently studied by several authors [124, 125]. Observational 
constraints have been determined by Sefusatti et al. [126]. Byrnes et al. performed a similar 
analysis in the special case of multiple-field models producing a local bispectrum [127, 128]. 
They allowed for deformations of the momentum triangle including a change of shape, but 
found these to be less important than rescalings of kf. 

Shape dependence is often substantially more complicated than scale dependence. Eq. (3.9) 
makes /nl dimensionless, but contains both powers and logarithms of the fcj. The powers 
occur as dimensionless ratios in which kt divides out, but the shape dependence remains. 
The argument of each logarithm is also a dimensionless ratio, but an extra scale is available: 



^^Chen implicitly worked in the equilateral limit ki — k, where kt — 3k and dlnfct — dlnfc. We are defining 
'^/nl tie the variation of /nl with perimeter for an arbitrary triangle if the shape is kept fixed. 
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the reference scale k-^. When present, this gives rise to the scahng logarithms In/sj/Zc^ and 
In kt/ki, described above, which depend on kt as well as the shape. It follows that a simple 
way to track the fet-dependence of /nl is to study the /c*-logarithms, yielding the identity 

For a general model, nj^^ can be written 

?^/nl = (guH^{hu - - v^,}fi{ki) + g2*{h2i. - v^}f2{ki) 



^{hs, -V,- 2s,}Mki) + ^{hi. - v,}h{ki) (3.30) 



where the fi{ki) functions are dimensionless ratios of polynomials in the ki which are listed 
in Table 14. 

Squeezed and equilateral limits. — In the equilateral limit, we find 



(3.31) 



The squeezed limit gives a simple result. 



n/NL (^g2^.{h2* - w*} + ^^{^3* - - v^}^ . (3.32) 



We define the running of the spectral index, Ug, by [122] 



Compatibility with (3.18) in the squeezed limit requires 

5 

In §5 we will verify this relation in the special case of canonical single-field inflation. 



4 Tensor modes 



Inflation will inevitably produce tensor fluctuations to accompany the scalar fluctuation 
Detection of the B-mode polarization signal produced by these fluctuations is a major aim 
of the Planck satellite and future CMB experiments. If it can be measured, this signal will 
provide important constraints on the energy scale of inflation. 

In certain models the tensor sector provides sufficient observables to allow one or more 
quantities, such as /nli to be written in terms of other observables. In the inflationary 
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literature such relationships are typically known as consistency relations, and were introduced 
by Copeland et al. [129, 130]. In the language of particle physics they are "observables in 
terms of observables" — predictions which are independent of how we parametrize the theory, 
and which the renormalization programme has taught us represent the physical content of 
any quantum field theory. Such consistency relations represent important tests of entire 
classes of models. To be used effectively with the next-order results of this paper we will 
require next-order predictions for the tensor modes. These were obtained by Stewart & Lyth 
[55], and are unchanged by the noncanonical action (1.1). 

In this section our aim is to obtain the next-order consistency relation, Eq. (4.6). The 
tensor fluctuation is a propagating spin-2 mode which belongs to the ADM field hij of (2.4). 



We write hij = a^e^^(e'^)jj, where ti jij 

1 
8 



0. At quadratic order, the action is [131] 



So 



d^x dr 



lijlij 



(4.1) 



There are two polarizations, traditionally denoted '+' and 'x,' making transverse in the 
sense di'yij = 0. Introducing a reference scale k^, and adding the power in each polarization 
incoherently, the resulting dimensionless spectrum can be written [55] 



2Hl 



1 + 2e J 1 - 7E - In 



2k 
kl 



(4.2) 



This is the sum of two copies of the power spectrum for a massless scalar field with = 1, and 
is conserved on superhorizon scales. Including next-order corrections, the scale dependence 
of Vg is measured by the tilt nt, 



dlnPg 
dln/c 



-2e. 



1 + - 7/^ 1 - 7E 



(4.3) 



It is conventional to measure the amplitude of tensor fluctuations relative to C,. One defines 
the tensor-to-scalar ratio, r, by the rule [130] 



^9 

where V is the dimensionless version of Eq. (2.22). We flnd 



16ej,c.o 



1 - 2??^ + (s^ + ??^) ( 7E + In 



2k 



(4.4) 



(4.5) 



In canonical models, r can be written purely in terms of observable quantities. In the 
noncanonical case this is not automatically possible without the addition of new observables. 
In general, 

( 1 2A;V 
, - r?^ + Sv, 7E + In — 



1 



k.J 



(4.6) 



One may use the lowest-order result for nt to eliminate e. To eliminate would require the 
scalar spectral index, n^. It is possible to use /nl to rewrite Cg in the prefactor [65], but 
in doing so one introduces dependence on the parameter I. Therefore at least two extra 
observables would be required to eliminate the dependence on s and If these depend on 
t, ^ or similar parameters, then further observables could be required. We conclude that at 
next-order, for a general P(X, (/>) Lagrangian, the observables {r, ?7-<j, nt, /nl} do not form a 
closed set. 
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5 Canonical single-field infiation 



The simplest model of inflation comprises a single scalar field with canonical kinetic terms. 
Maldacena showed that the fluctuations in this model are almost Gaussian, with /nl of order 
r [26]. This is unobservably small. In a canonical model, Cg = 1. 

Nonlinearity parameter. — To calculate /nl we require the flow parameters hi, which measure 
time dependence in the vertex factors gi. These are 

T]{2e -r]-0 

] (5.1) 



hi=0 h2 



e — ri 



/i3 = : /i4 = z h = 3ri- ^ 

e + rj 4 

In this model the time-dependence of z is described hy v = rj. Collecting contributions from 
Table 5, the equilateral limit of /nl can be written 



f 5 
Jnl 1^ 

00 



11^^ + + ^ (768^ - 54) + I37E - 8 + 3 In 

' 216 \ J 36 \ ' K 

+ ^%^|ll7E-14 + 64a. + llln^ 
36 K 



(5.2) 



where we have used the numerical constant a; = coth^^ 5, and k should be regarded as the 
common momentum scale, ki = k. The squeezed limit may be recovered from Table 6. We 
find 



/nl^ Y2 



Ik) ( 2k)' 

2e^ + ri^ + 2ei + r/^C^<j 7e - 2 + In — ^ + e^r]A 27E - 1 + 21n — ^ 



(5.3) 



where k should now be regarded as the scale of the hard momenta in the correlation func- 
tion. In §3.4 and §3.6 we emphasized that /nl should be finite in this limit, containing no 
large logarithms, because these factorize into the power spectrum and are subtracted. The 
remaining logarithms [the ln2/c/A;* terms in (5.3)] track the dependence of /nl on the hard 
scale, and will be studied below. Using (3.18) and comparing with the spectral indices quoted 
in Table 1, it is easy to check that our formula correctly reproduces the Maldacena limit. 
We note that, strictly, one should regard agreement in this limit as an accident which hap- 
pens because the simple slow-roll model contains no other scale which could interfere with 
factorization of the correlation function. 

Scale dependence of f^i,. — Specializing to the equilateral limit of n^^^, we find 

^/nl ^ ^^*(66e- + 18^*). (5.4) 
In the squeezed limit one obtains 

5 5 

which correctly describes the running of the scalar spectral index, a<j*, in agreement with (3.34). 
The consistency conditions (5.3) and (5.5) represent a nontrivial check on the correctness of 
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our calculation. In particular, throughout the calculation we have cleanly separated the 
conceptually different scales kt and k^. Therefore the correct formula (5.5) is not simply 
a consequence of obtaining the correct lowest-order terms in (5.3) — although, as described 
in §2.2, it can be obtained from these by differentiating with respect to In/c^ after setting 
ki, = kt. 



6 Non-canonical single-field infiation 

For the noncanonical action (1.1), the hi can be written 



hi=e + r]-2s + — — ^ — ^ = r/ - 4s + 



iqi_2A ^ ' ■^_e-3(l-c2) 



r?(6 + g)-2s(t-c^ ) , Ve \ ^^-l) 

e + ?? - 2s + (1 - 



/i3 = r/ - 2s + ; K /i4 = 2r? - 4s 

e + ?? — 2s + (1 — H 

h5 = 3r]- 4s, 



where we have defined ^ = r]/Hrj, t = s/Hs. 

In the canonical case, it was possible to verify Maldacena's consistency condition to next- 
order. In the noncanonical case this is not possible without a next-next-order calculation, 
because for 7^ 1 the leading contribution to /nl is 0(1) in the slow-variation expansion. 
Therefore our calculation of subleading corrections produces a result valid to 0(e), which is 
short of the O(e^) accuracy required to verify the consistency condition at next-order. Chen 
et al. gave the subleading corrections in terms of undetermined integrals [43]. Expanding 
these asymptotically, they argued that the consistency relation would be satisfied at lowest- 
order. More recently, Renaux-Petel [114] gave an equivalent demonstration. Here, we have 
knowledge of the full bispectrum to subleading order. Using (6.1), it can be verified that in 
the squeezed limit, and expanding around a reference scale A:*, 

/nl ^ ^ (2e^ + r?^ + s^ . (6.2) 

One may check that this agrees with Eq. (3.18) and Table 1. We expect n/j^^ = O(e^), 
and therefore a next-next-order calculation is required to estimate the running of /nl in 
noncanonical models. 



6.1 Asymptotically power-law models 

Power-law inflationary models were introduced by Lucchin & Matarrese [132, 133], who 
studied potentials producing an expansion history of the form a{t) oc t^/^ . The exponent 1/e 
is the usual parameter e = —H/H"^. It need not be small, but should be taken as constant 
which makes = ^ = Q. The solution is inflating provided e < 1. Exact solutions can be 
found in the canonical case, which form the basis of the next-order calculation [60]. 

In this section we study two examples which are asymptotically described by noncanon- 
ical power- law inflation at late times. The flrst is Dirac-Born-Infeld ( "DBI" ) inflation, which 
produces a scale-invariant power spectrum at lowest-order. Departures from scale invariance 
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appear at next-order. These properties imply that we can compare our results to a formula 
of Khoury & Piazza which was obtained without invoking the slow-roll approximation [134]. 
Our second example is ^-inflation, for which the power spectrum is not scale invariant at 
lowest-order, and to which Khoury & Piazza's result does not apply. 

6.1.1 Dirac Born Infeld inflation 

The DBI action is a low-energy effective theory which describes a D3-brane moving in a 
warped throat. It was proposed as a model of inflation by Silverstein & Tong [61], and 
subsequently developed by Alishahiha et al. [33]. The action is of the form (1.1), with 
P{X, (/)) satisfying 

PiX, ^) = -j^^[y/l- f{ct>)X - l] - V{(p) (6.3) 

where / is an arbitrary function of (p known as the warp factor, and V{(t)) is a potential arising 
from couplings between the brane and other degrees of freedom. The DBI Lagrangian is 
algebraically special [65, 125, 135] and enjoys a number of remarkable properties, including a 
form of nonrenormalization theorem [136-138]. In principle non- minimal curvature couplings 
can be present, of the form RcjP', which spoil inflation [139]. This gives a form of the rj- 
problem, and we assume such terms to be negligible. 

Eq. (6.3) makes 2A/S = (1 - cl)/cl, which requires g'l — )• but causes the denominator 
of hi in (6.1) to diverge. Only the finite combination gihi appears in physical quantities, and 
it can be checked that gihi — )• as required. The square root in (6.3) must be real, giving a 
dynamical speed limit for (j). It is conventional to define a Lorentz factor 

7 = (l-/<^')-^/'. (6.4) 

When 7 ~ 1 the motion is nonrelativistic. When 7^1, the brane is moving close to the 
speed limit. The Lorentz factor is related to the speed of sound by Cs = 7"^. 

Silverstein & Tong [61] argued that (6.3) supported attractor solutions described at 
late times by power-law inflation. In this limit, the slow-variation parameters e and s are 
constant, with 7/ = ^ = t = but £ not zero. Variation of the sound speed gives s = — 2e, 
making = and yielding scale-invariant fluctuations at lowest-order [cf. (2.23)]. In the 
equilateral limit, ^'^ we find that /nl satisfies 

v2 



/nl ~ 1) 



l-^^(3-47E)e + 0(7.-^) 



(6.5) 



In §1 we estimated the relative uncertainty in /nl to be ~ 14e, working in the limit 7^1, 
based on 0(e) terms from the vertices only. Eq. (6.5) shows that, due to an apparently 
fortuitous cancellation, this large contribution is almost completely subtracted to leave a 
small fractional correction ~ 0.69e. 
In the squeezed limit we find 

/NL^^e2(47E-5 + 41n|). (6.6) 



^^Recall that the equilateral limit is /NL(fc,fc,fc) and is not the quantity constrained by experiment. See 
footnote 17. 
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which is O(e^), as predicted by the Maldacena condition (3.18) and the property w = [33]. 

Comparison with previous results. — Khoury & Piazza estimated the bispectrum in a power- 
law model satisfying zu = without invoking an expansion in slow-variation parameters 
[134].^^ They quoted their results in terms of a quantity fx which replaces A, 



For the DBI model, fx = 1 — c^. They assumed constant fx, making their result valid to all 
orders in e but only lowest-order in the time dependence of fx- Working in the equilateral 
limit for arbitrary constant fx we find 



/nl^-7t;^[55(1-cL) + 8/x] 



972c2^ 

+ 



972c2 



s-k 



3k I 3k 

149 - 8cl^ - 2207E - 220 In — + /x <^ 40 - 327e - 32 In — 



(6.8) 



Adopting the evaluation point k^ = 3k, this precisely reproduces (8.4) of Khoury & Piazza 
[134] when expanded to order e. Although (6.8) does not strictly apply to DBI, where Cg is 
changing, it can be checked that effects due to the time dependence of fx do not appear at 
next-order. Indeed, Eq. (6.8) yields (6.5) when /x = 1 — 7~^- 

Generalized DBI inflation. — The foregoing analysis was restricted to the asymptotic power- 
law regime, but this is not required. Using an arbitrary potential V{(j)) in (6.3) one can 
construct a generic quasi-de Sitter background. Many of their properties, including the 
attractor behaviour, were studied by Tranche et al. [141]. However, in the absence of a 
controlled calculation of next-order terms it has not been possible to estimate corrections from 
the shape of V{4>) or /(</>). Analogous effects have been computed for Galileon inflation [34], 
but our computation enables us to determine them in the DBI scenario for the first time. 
For 7^1 the noncanonical structure suppresses background dependence on details of the 
potential. But small fluctuations around the background cannot be shielded from these 
details, which induce three-body interactions whether or not they are relevant for supporting 
the quasi-de Sitter epoch. These interactions generate relatively unsuppressed contributions 
to the three-point function. 

We adapt the notation of Franche et al., who defined quantities measuring the shape of 
the potential V and the warp factor /,22 



2Vy;' --y ^ = ^sn(0^')7fe3^- (^-^^ 

The same branch of /^/^ should be chosen in computing /'^/^ ^nd sgn{4>f^^'^). Note that 
these shape parameters do not coincide with the global slow- variation parameters e and r/. 



^^See also Baumann et al. [140]. 

^^The factor sgn(</>/^''^) was not used by Franche et al. [141], but is necessary here because the relativistic 
background solution requires fX = 1 + 0(7~'^) up to corrections suppressed by 0(e) which are higher-order 
than those we retain. Depending on the direction of motion, this yields = ±/~^^'^ + 0(7~^). 
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Franche et al. argued that A <C 1 was required to obtain attractor solutions, which we will 
assume to be satisfied in what follows. In addition, we work in the equilateral limit and take 
7 S> 1, which is the regime of principal interest for observably large /nl- We find 

£ = — , 7] = 3 A, and s = ! A. (6.10) 

7 772 772 

The leading term of (6.5) is unchanged, but the subleading terms are dominated by shape- 
dependent corrections. 



, 357 

JNL = 



108 



3A,^ / , 2k\ ri^^ ( 3k\ 

1 + — ^ 31 + 147E- 228a; + 141n— + ;P 3 - 147e - 141n — 
14 V / 77^ V fc^ / 

U3 - 77E - 256a; - 71n ^) + 0(7"') 



77v, V k^J 



(6.11) 



These subleading terms are more important than those of (6.5), which began at relative order 
7~-^ and are therefore strongly suppressed for 7 S> 1. Moreover, inflation can occur even for 
relatively large values of and r]y — roughly, whenever ev/7 < 1 — so these corrections need 
not be extremely small. For large |/nl|) we estimate the relative correction to be 

A/nl_ , 2.10g.. - 0.41r?., 

— ss-2.75A,, + . (6.12) 

JNL I/nlT/^ 

For negative equilateral /nl, current constraints approximately require |/nl|^''^ ^ 12 [6]. 
Therefore, these corrections can be rather important unless the potential is tuned to be flat, 
although some cancellation occurs because and r/v enter with opposite signs. 

To obtain an estimate, suppose that A^ is negligible. Taking the extreme 95%-confidence 
value /nl = — 151 [6] and ~ |^v| ~ 1 to obtain an estimate for a "generic" potential, the 
correction is of order 14% if 77^ > and 20% if ??v < 0. To reduce these shifts one might 
be prepared to tolerate a small tuning, giving perhaps ~ \Vv\ ~ 0.1 and suppressing 
the correction to the percent level. However, the corrections grow with decreasing |/nl|- 
Keeping the generic estimate ~ \r]v\ ~ 1, and using |/nl| ~ 50, for which 7 10 and 
the approximation 7^1 used to derive (6.12) is at the limit of its applicability, we find the 
corrections to be of order 24% for r]^ > and 36% for r/v < 0. 

Although Franche et al. argued that A^ must be small to obtain attractor behaviour, it 
need not be entirely negligible. In such cases it introduces a dependence on the shape of the 
warp factor in addition to the shape of the potential. This may be positive or negative. If 
the A^ and terms add constructively, the next-order correction can become rather large. 

Infrared model. — The DBI scenario can be realized in several ways. The original "ultraviolet" 
model is now disfavoured by microscopic considerations [63, 65]. Chen introduced [123, 142, 
143] an alternative "infrared" implementation which evades these constraints and remains 
compatible with observation [66, 144-146]. In this model the warp factor /(</>) is A/(/)^, 
where A is a dimensionless parameter. The potential is V{(p) = Vq — fS'^H'^cp'^ /2, in which the 
mass is expressed as a fraction /3^/^ of the Hubble scale. The constant term Vq is taken to 
dominate, making e negligible. However the remaining slow-variation parameters need not 
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be small. It is convenient to express our results in terms of the number, Nf,, of e- folds to the 
end of inflation. Background quantities evaluated at this time are denoted by a subscript 
'e'. Computing (6.10) with these choices of V and / we find rye ^ S/A'^g and Sg ~ l/Ng. 
The infrared model is an example where Ag is not negligible, being also of order 1/Ne- 
Specializing (6.12) to this case, we find 

^.-^(o..H..-.M..:.>n|).^, ,e.:3) 

where we have chosen k-^ = 3k in the final step. Adopting the best-fit value A'^e ~ 38 
suggested by Bean et al. [66], we find a fractional correction of order 12%. This relatively 
small correction is a consequence of the negligible e in this model. The analysis of Bean 
et al. gave a reasonable fit for a range of /3 of order unity. Keeping ~ 38 and using 
the maximum likelihood value /3 = 1.77 quoted by Bean et al., we find A/nl ~ —19. The 
corresponding shift is from /nl ~ —163 without next-order corrections to /nl ~ —182 with 
next-order corrections included. 

6.1.2 A;-infiation 

The A;-inflation model of Armendariz-Picon et al. [39] also admits power-law solutions. The 
action is 

4 4 - 37 - X 
9^ 

where 7 is a constant, no longer related to the speed of sound by the formula Cs = 7~^ which 
applied for DBI. Unlike the DBI Lagrangian, Eq. (6.14) is unlikely to be radiatively stable 
and its microscopic motivation is uncertain. Nevertheless, nongaussian properties of the 
inflationary fluctuations in this model were studied by Chen et al. [43]. There is a solution 
with 

making e = 3^/2 and Cg constant. Therefore this model has s = but e 7^ 0, and is not 
scale-invariant even at lowest-order. Inflation occurs if < 7 < 2/3. The lowest-order 
contribution to /nl is of order I/7, making the next-order term of order unity. A next-next- 
order calculation would be required to accurately estimate the term of order 7. 

In the equilateral limit, Chen et al. quoted the lowest-order result /nl ~ — I7O/8I7. 
Proceeding as in §1, one can estimate the fractional theoretical uncertainty in this prediction 
to be ~ 97, or roughly ±20. This is comparable to the Planck error bar, and is likely to 
exceed the error bar achieved by a subsequent CMB satellite. Still working in the equilateral 
limit, we find 

170 



PiX,<p) = - / , (6.14) 



8I7 



1-^(61-I921n0+O(72) 



(6.16) 



As for DBI inflation, a fortuitous cancellation brings the fractional correction down from our 
estimate ~ 97 to ~ O.57. It was not necessary to choose a reference scale /c* in order to 
evaluate (6.16). In DBI inflation, to the accuracy of our calculation, the power spectrum 
is scale invariant but /nl is not. For the power-law solution of fc-inflation, with the same 
proviso, it transpires that /nl is scale invariant even though the the power spectrum is not. 
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Because w ^ in this model, the analysis of Khoury & Piazza does not apply. In a 
recent preprint, Noller & Magueijo gave a generalization which was intended to be valid for 
small w and constant but otherwise arbitrary e and s [147]. Their analysis also assumes 
constant fx- We set the reference scale to be k^, = 3k and work in the equilateral limit for 
arbitrary constant fx- One finds 

/nl^-^^[55(1-cL) + 8/x] 

+ 9^ [177 + 120cL - 10240.(1 - 4) + /x{264 - 1280^}] (6.17) 

+ [7 + 557E + 324 - 256^^(1 - 4) + /x{56 + 87E - 320a;}] . 

For s = — 2e, both (6.17) and Noller &: Magueijo's formula (A. 15) reduce to (6.8) evaluated 
at k-^ = 3k. For s 7^ — 2e, Eq. (6.17) disagrees with Noller &; Magueijo's result. This occurs 
partially because they approximate the propagator (2.20) in the superhorizon limit |A;Cst| ^ 1 
where details of the interference between growing and decaying modes around the time of 
horizon exit are lost. For example, their approximation discards the Ei-contributions of (A. 2) 
although these are 0{zu) and as large as other contributions which are retained. But were 
these terms kept, the superhorizon limit |A;Cst| ^ 1 could not be used to estimate them. 
Infrared safety of the Ji integrals in (3.15) is spoiled if truncated at any finite order, causing 
divergences in the squeezed limit iJj — t- and a spurious contribution to the bispectrum 
with a local shape. As we explain in Appendices A and B, it appears that — as a point of 
principle — if tz7 / corrections are kept then the shape of the bispectrum can be accurately 
determined only if the full time-dependence of each wavefunction around the time of horizon 
exit is retained. 



7 Conclusions 

In the near future, we can expect key cosmological observables to be determined to high 
precision. For example, the Planck satellite may determine the scalar spectral index Ug to an 
accuracy of roughly one part in 10^ [148]. In a formerly data-starved science, such precision 
is startling. But it cannot be exploited effectively unless our theoretical predictions keep 
pace. 

Almost twenty years ago, Stewart & Lyth developed analytic formulae for the two- 
point function accurate to next-order in the small-parameter e = —H/H^. Subsequent 
observational developments have restricted attention to a region of parameter space where 
e <C 1 is a good approximation, making the lowest-order prediction for the power spectrum 
an accurate match for experiment. The same need not be true for three- and higher n-point 
correlations, where the imminent arrival of data is expected to improve the observational 
situation. The results of §6 show that Planck^s observational precision in the equilateral mode 
may be comparable to next-order corrections. For a future CMB satellite it is conceivable 
that the data will be more precise than a lowest-order estimate. In this paper we have 
reported a next-order calculation of the bispectrum in a fairly general class of single-field 
inflationary models: those which can be described by a Lagrangian of the form P{X,(j)), 
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where X = g 'Va4>'^b4' is twice the field's kinetic energy. Our calculation can be translated 
to /nL) and in many models it provides a much more precise estimate than the lowest-order 
result. 

Our major results can be categorized into three groups. The first group comprises tests 
of the accuracy, or improvements in the precision, of existing lowest-order calculations. The 
second generate new shapes for the bispectrum at next-order. The third involves technical 
refinements in calculating scalar n-point functions for n ^ 3. 

7.1 Accuracy and precision 

Except in special cases where exact results are possible — such as the result of Khoury & Piazza 
for constant fx discussed in §6 — predictions for observable quantities come with a "theory 
error" encapsulating uncertainty due to small contributions which were not calculated. For 
inflationary observables the typical scale of the theory error is set by the accuracy of the 
slow-variation approximation, where the dimensionless quantities e = —H/H^, t] = e/He, 
s = Cs/Hcs (and others) are taken to be small. 

Power-law DBI inflation and k-inflation. — In §1 and §§6.1.1-6.1.2 we estimated the preci- 
sion which could be ascribed to the lowest-order formula for /nl in the absence of a com- 
plete calculation of next-order effects. To do so one may use any convenient — but hopefully 
representative — subset of next-order terms, estimating the remainder to be of comparable 
magnitude but uncertain sign. Using the next-order contributions from the vertex coef- 
ficients Qi, which can be obtained without detailed calculation, we estimated the fractional 
uncertainty to be of order 14e for DBI and 97 for /c-inflation. The next-order terms neglected 
in this estimate come from corrections to the propagator, and from the time-dependence of 
each vertex. The prospect of such large uncertainties implies one has no option but to carry 
out the full computation of all next-order terms. 

In the power-law DBI and A;-inflation scenarios, we find that the terms omitted from 
these estimates generate large cancellations, in each case reducing the next-order contribution 
by roughly 95%. Similar large cancellations were observed by Gong & Stewart in their 
calculation of next-next-order corrections to the power spectrum [57]. After the fact, it 
seems reasonable to infer that the contributions from gi systematically overpredict the next- 
order terms. But this could not have been deduced without a calculation of all next-order 
effects. Therefore, for power-law DBI and /c-inflation models we conclude that the lowest- 
order calculation is surprisingly accurate. 

For DBI inflation the status of the power-law solution is unclear, being under pressure 
from both observational and theoretical considerations. More interest is attached to the 
generalized case to be discussed presently. For A;-inflation, taking present-day constraints on 
the spectral index into account, the next-order correction is of order 1%. Estimating the 
contribution of next-next-order terms using all available contributions from our calculation 
gives a fractional uncertainty — measured with respect to the lowest-order term — of order 
227^^. If similar large cancellations occur with terms not included in this estimate, the next- 
order result could be rather more precise than this would suggest.^"' Without knowledge of 

^^As suggested in the introduction, these terms could be calculated using the next-next-order propagator 
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such cancellations, however, we conclude that the uncertainty in the prediction for /nl has 
diminished to ~ 2507% ^ W% of the uncertainty before calculating next-order terms. In 
the case of power-law DBI inflation, the same method yields an estimate of next-next-order 
corrections at ~ 40e^ , also measured from the dominant lowest-order term in the limit 7^1. 
This reduces the uncertainty to ~ 300e% < 15% of its prior value. 

Generalized DBI inflation. — The situation is different for a generalized DBI model with ar- 
bitrary potential V{(j)) and warp factor /(</>). The largest next-order corrections measure a 
qualitatively new effect, not included in the power-law solution, arising from the shape of V 
and /. The fractional shift was quoted in (6.12) and can be large, because the DBI action 
supports inflation on relatively steep potentials. Indeed, if one were to tune the potential 
to be flat in the sense Ev ~ |f?v| ^ 10~^ then much of the motivation for a higher-derivative 
model would have been lost. Even for rather large values of |/nl| the correction can be sev- 
eral tens of percent for an "untuned" potential with Ev ~ |^v| ~ 1- For slightly smaller |/nl| 
the correction is increasingly significant, perhaps growing to ~ 35%. The formulas quoted 
in §6.1.1 assume 7 S> 1 and would require modification for very small /nl where 0(7"^) 
corrections need not be negligible. If desired, these can be obtained from our full formulae 
tabulated in §3.2. 

In a concrete model — the infrared DBI scenario proposed by Chen [123, 142, 143] — we 
find the correction to be ~ 12% for parameter values currently favoured by observation, which 
translates to reasonably large shifts in /nl- For the maximum-likelihood mass suggested 
by the analysis of Bean et al., we find next-order corrections increase the magnitude of 
/nl by a shift |A/nl| ~ 19. This is a little smaller than the error bar which Planck is 
expected to achieve, but nevertheless of comparable magnitude. We conclude that a next- 
order calculation will be adequate for Planck, but if the model is not subsequently ruled out 
a next-next-order calculation may be desirable for a CMBPol- or CoRE-type satellite. 

For models producing small /nlj such as the "powerlike" Lagrangian discussed in §VI.B 
of Franche et al. [141], we find similar conclusions. However, in such models /nl is unlikely 
to be observable and the subleading corrections are of less interest. 

7.2 New bispectrum shapes 

Because our final bispectra capture the shape dependence in the squeezed limit, we are able 
to determine the relationship between the lowest-order and next-order shapes, discussed in 
§3.5. 

Working in a P{X, (p) model, the enhanced part of the lowest-order bispectrum is well- 
known to correlate with the equilateral template. Only two shapes are available, plotted in 
Table 9, and the bispectrum is a linear combination of these. The next-order bispectrum is a 
linear combination of seven different shapes, although these cannot be varied independently: 
strong correlations among their coefficients are imposed by the P{X, <j)) Lagrangian. Many 
of these shapes also correlate with the equilateral mode, but two of them are different: in 
the language of §3.5, these are S'^ and Ss- For typical values of a (which is the ^-enhanced 

corrections provided by Gong & Stewart [57]. 
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part of A/S) one can obtain a linear combination of the next-order shapes orthogonal to both 
lowest-order shapes. This is the orthogonal shape O appearing in Table 11. 

This shape represents a distinctive prediction of the next-order theory, which cannot be 
reproduced at lowest-order. It is very similar to a new shape constructed by Creminelli et 
al. [37], working in an entirely different model — a Galilean-shift invariant Lagrangian with 
at least two derivatives applied to each field. We conclude that even a clear detection of a 
bispectrum with this shape will not be sufficient, on its own, to distinguish between models 
of P{X, (p)- and Galileon-type. In practice, a network of interlocking predictions is likely to 
be required. For example, this shape cannot be dominant in a P{X, (/>) model. It must be 
accompanied by a mixture of Si- and 5'2-shapes of larger amplitude. 

The orthogonal shape is not guaranteed to be present at an observable level in every 
model, even in models where the next-order corrections are detectable. This may happen 
because a takes a value for which no bispectrum orthogonal to Si and 52 can be generated 
perturbatively, or because the slow-variation parameters accidentally conspire to suppress its 
amplitude. Nevertheless, it will be present in many models. If detected, it would play a role 
similar to the well-known consistency relation between r and ut. Searching for the 0-shape 
in real data is likely to require a dedicated template, and until this is constructed it is not 
possible to estimate the signal-to-noise and therefore the amplitude detectable by Planck or 
a subsequent experiment. It would be interesting to determine the precise amplitude to be 
expected in motivated models, but we have not attempted such an estimate in this paper. 

7.3 Technical results 

Our calculation includes a number of more technical results. 

Treatment of boundary terms. — In §3.1 we gave a systematic treatment of boundary terms 
in the third-order action. Although these terms were properly accounted for in previous 
results [26, 42, 43], these calculations used a field redefinition which was not guaranteed to 
remove all terms in the boundary action. 

Pure shape logarithms. — The subleading correction to the propagator contains an exponential 
integral contribution [of the form Ei(x)] whose time dependence cannot be described by 
elementary functions [34, 43, 57]. This term must be handled carefully to avoid unphysical 
infrared divergences in the squeezed limit where one momentum goes to zero. In Appendix B 
we describe how this contribution yields the Jj functions given in Eq. (3.15). These are 
obtained using a resummation and analytic continuation technique introduced in Ref. [34]. 
The possibility of spurious divergences in the squeezed limit shows that, as a matter of 
principle, one should be cautious when determining the shape of the bispectrum generated 
by an approximation to the elementary wavefunctions. 

In the present case, the Jj contain 'pure' shape logarithms which are important in 
describing factorization of the three-point function in the limit fcj — t- 0. Obtaining the 
quantitatively correct momentum behaviour requires all details of the interference in time 
between growing and decaying modes near horizon exit. The possibility of such interference 
effects, absent in classical mechanics, is a typical feature of quantum mechanical processes. 
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This interference correctly resolves the unwanted divergences in the opposing infrared limit 
kj — )• with j ^ i- We discuss these issues in more detail in Appendices A and B. 

Comparison with known results. — We have verified Maldacena's consistency condition (3.18) 
to next-order in canonical models by explicit calculation of the full bispectrum. This agrees 
with a recent calculation by Renaux-Petel [114]. In the case of power-law inflation with 
vj = and constant e and s, we reproduce a known result due to Khoury & Piazza [134]. 

A subset of our corrections were computed by Chen et al. [43]. Up to first-order terms in 
powers of slow-roll quantities — where our results can be compared — we find exact agreement. 

Note added 

Immediately prior to completion of this paper, a preprint by Arroja & Tanaka appeared 
[149] which appears to present arguments regarding the role of boundary terms which are 
equivalent to those of §3.1.2. 
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operator 










/nlIio 


5 9i* 


5 92* 
8 z,c2. 




l^h\i 


, 2k 

7e + In — 

0.577216'^ 


4 , 2k 
7E-3+ln- 

-0.756118" 






2k 

-7e + 1 ^ In 
0.422784'^ 


5 , 2fc 
1.08945'' 






3 


2.51224^^ 






2 


2 

3 

0.666667 





° Evaluated at the reference scale fc^ = 2A:uvj where fcuv is the common hard momentum 



Table 6. Squeezed limit of /nl at lowest-order and next-order. The numerical constant uj satisfies 
w = 5 In I = coth^^ 5. 













Ss 


S's 


S'g 


locaP 




0.50*= 


0.37" 


OAS" 


0.54"= 


0.39" 


0.42*= 


equilateral^ 


0.99 


0.87 


1.00 


0.93 


0.80 


0.94 


0.94 


orthogonaF 


0.084 


0.46 


0.065 


0.31 


0.52 


0.25 


0.29 


enfolded'^ 


0.60 


0.86 


0.59 


0.77 


0.87 


0.72 


0.75 



Sec Komatsu Sz Spergel [8] and Babich ct al. [118]. 
Sec Babich ct al. [118]. 

^ See Senatorc ct al. [108]. 

See Mecrburg ct al. [121] and Senatorc ct al. [108]. 

^ The local template, and the operators CC'^ and C(^C)^! strongly peaked in the "squeezed" 
limit where one momentum becomes much softer than the other two. For these shapes the 
inner product which defines the cosine is divergent, and must be regulated. The resulting 
cosines are almost entirely regulator-dependent. See the discussion in §3.5. 

The values we quote are meaningful only for our choice of regulator. For the values quoted 
above we have used (5j^^j-^ — k/kf — 10~^, where <5min was defined in the main text. 



Table 7. Overlap cosines for the bispectrum shape proportional to each slow-variation parameter. 
Sign information has been discarded. 
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local 


equilateral 


orthogonal 


enfolded 


local 


1.00 








equilateral 


0.34 


1.00 






orthogonal 


0.49 


0.03 


1.00 




enfolded 


0.60 


0.51 


0.85 


1.00 



Table 8. Overlap cosines between common templates, defined in Table 7. 



Babich" et al. Fergusson & Shellard'' 



Shape 1 



Shape 2 



^ Sec Babich ct al. [118]. The plotted quantity is x'^ y"^ B(x , y , 1) , where x — ki/k^, y = 
^2/^3 and B is the bispectrum, and normalized to unity at the equilateral point x = y — 1. 

^ See Fergusson & Shellard [119]. The plotted quantity is fc^fc2^|B(fci, k2, ^3) as a function 
of the CK and j3 parameters defined in {3.23a)— {3.23c) . 



Table 9. Lowest-order bispectrum shapes enhanced by ^ in P{X,(f)) models. 
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Babich et al. 



Fergusson & Shellard 




Table 10. Bispectrum shapes enhanced by ^ at next-order in a P{X, (p) model. The Babich et al. 
and Fergusson-Shellard plots are defined in Table 9. 
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Babich ct al. 



Fergusson & Shellard 



a = 10-3 




Table 11. The orthogonal bispectrum shape O has zero overlap with both the lowest-order possi- 
bilities S'l and 5*2. The Babich et al. and Fergusson-Shellard plots are defined in Table 9. 



Babich ct al. Fergusson & Shellard 




Table 12. Highly orthogonal shape constructed by Creniinelli et al. [37]. The Babich et al. and 
Fergusson-Shellard plots were defined in Table 9. 
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a 


local 


equilateral 


orthogonal 


enfolded 


Creminelli et al.'' 




0.35" 


0.012 


0.36 


0.32 


0.89 


1 


o.as'^ 


0.012 


0.38 


0.33 


0.86 


10 


0.4P 


0.011 


0.40 


0.35 


0.81 



^ The local template is divergent, and the values we quote arc meaningful only for our choice of 
regulator. For the values quoted above wc have used (ij-^in — k/kt — 10~'^, where '5n-^il-^ was 
defined in §3.5. 

^ For the Creminelli et al. shape, see Table 12. 



Table 13. Overlap cosines between the orthogonal shape O and common templates, defined in 
Table 7. Sign information has been discarded. 
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A Propagator corrections 



At lowest-order, the wavefunctions can be obtained from Eq. (2.20) by setting all slow- 
variation parameters to zero. Choosing a reference scale A;^, this yields the standard result 

Uk{r) = i \ (1 - ikcs.r)e''^-^. (A.l) 

Next-order corrections to the propagator were discussed by Stewart & Lyth [55], who 
quoted their result in terms of special functions and expanded uniformly to next-order after 
taking a late-time limit. The uniform next-order expansion at a generic time was given by 
Gong & Stewart [57] for canonical models and by Chen et al. [43] in the noncanonical case. 
Their result was cast in a more convenient form in Ref. [34] whose argument we briefly review. 

The next-order correction is obtained after systematic expansion of each quantity in 
(2.20) to linear order in the slow-variation parameters. Contributions arise from each time- 
dependent factor and from the order of the Hankel function. Collecting the formulae quoted 
in Ref. [43], one finds 



dH^ix) 



dv 



i /2 



3.3/2 y ^ 

^=2 



e'^(l-ix)Ei(-2i2;)-2e-'^-i^e-'^(l + ix) , (A.2) 



where Ei(x) is the exponential integral, 

Ei(x) = / -di if X G M. (A.3) 

J —00 ^ 

This is well-defined for x < 0. For x ^ it should be understood as a Cauchy principal value. 
For complex argument — required by (A.2) — one promotes Eq. (A.3) to a contour integral, 
taken on a path running between t = x £ C and \t\ — )• 00 in the half-plane Re(t) < 0. Using 
Cauchy's theorem to rotate this contour onto the negative real t-axis, one finds 

Ei(2iA:c,r) = lim ^ e^'^^^^. (A.4) 

i-oo(l+ie) C 

The next-order correction to (A.l), expanded uniformly to 0(e) in slow-variation parameters 
but including the exact time-dependence, is 



5uik, t) = ^ * <^ - ti7.e-"=^-^(l + ikcs.T) / -f e 

/xo* + inukcsi^T + s^^k^cl^r"^ + AiV^ (tUj, - itx^^fcc^^r - Sy^/c^cf^r^) 



(A.5) 



where 



vr vr 
^LQ = e + V + 2s + \—w and ^i = e + s — i— ro, (A. 6) 

and quantities labelled are evaluated at the horizon-crossing time for the reference scale 
ki,. We have used AA^^ to denote the number of e- folds which have elapsed since this time. 
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so AA^* = In |/cjfCs*r|. The limit e — )• should be understood after the integration has been 
performed and merely guarantees convergence as ^ — )• — oo. 

Setting the integral aside, Eq. (A. 5) has the appearance of an expansion in powers of kr. 
In fact, the time dependence is exact. Although terms of high-order in kr become increasingly 
irrelevant as |A;r| — t- oo, all such terms make comparable contributions to a generic n-point 
function at the time of horizon crossing, where |A:r| ~ 1, and should not be discarded. This 
reflects interference effects between the growing and decaying modes at horizon exit. Indeed, a 
high-order term such as (/cr)" typically generates a contribution ~ n\{k/kt)^. Such terms are 
suppressed by the factor {k/kt)"^ which typically varies between and 1/2, but are enhanced 
by the rapidly growing factorial. Therefore their contribution to the shape dependence must 
usually be retained. An infinite series of such terms may converge only for certain ratios 
k/kt, requiring a continuation technique to obtain the momentum dependence — and hence 
the shape of the bispectrum — for arbitrary ki . We will see an explicit example in Appendix B 
below. 

To evaluate the cubic action (3.10) we will require the time derivative 5u' , where ' 
denotes a derivative with respect to conformal time. Using the identity ^lo + m — 2wi = 0, 
it follows that 



In order that the power spectrum remains conserved on superhorizon scales, one should 
verify that 5u approaches a constant and 5u' — )• as |A;r| — t- 0. The explicitly time-dependent 
term AA^* apparently spoils the required behaviour, but is compensated by a logarithmic 
divergence from the integral. Indeed, for r — t- one finds 



which precisely cancels the time dependence arising from AN^,. Note the incomplete can- 
cellation of the logarithm, which leaves a residual of the form zu^. In 2 and is the origin of 
the In 2 term in the Stewart-Lyth constant C = —2 -|- ln2 -|- 7e [55, 60]. This is a primitive 
form of the incomplete cancellation which leads to residual In ki / kt and In kt / fc* terms after 
cancellation of In r logarithms in the three-point function. 

B Integrals involving the exponential integral Ei(z) 

The principal obstruction to evaluation of the next-order corrections using standard methods 
is the Ei-term in (A. 5). This can not be expressed directly in terms of elementary functions 
whose integrals can be computed in closed form. It was explained in Appendix A that 
although Ei(2iA;Csr) contains terms of high orders in kr which become increasingly irrelevant 
as |/cr| —7- 0, these cannot usually be neglected when computing n-point functions. A term of 
order (A;r)" generates a contribution ~ n\{k/kt)'^ and as part of an infinite series will converge 




(A.7) 




(A.8) 
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at best slowly for certain values of k/kt — typically, when the scale k itself is squeezed to zero. 
The ratio k/kt is maximized when k is left fixed but one of the other momenta is squeezed. 
To obtain a quantitative description of this limit may require continuation. Hence, in order 
to understand the shape generated by arbitrary ki, it follows that the entire A;r-dependence 
must be retained. It was for this reason that we were careful to expand uniformly only in 
slow-roll parameters in Appendix A, while keeping all time-dependent terms. 
The integrals we require are of the form [34] 



We are again using the convention that, although depends on all three ki, only the asym- 
metric momentum is written as an explicit argument. In the following discussion we specialize 
to Iq. Comparable results for arbitrary m can be found by straightforward modifications and 
are given in Ref. [34]. 

We introduce the dimensionless combinations 

After contour rotation and some algebraic simplification, we can express Iq in terms of tJs 
and 6*3 [34] 

-u / -2v 



Ioik3) = -^T— / e-" / - e-'\ (B.3) 

VshCs Jo 



The v-integral has a Puiseaux series representation 

where the sum converges uniformly for all complex O^u. The right-hand side has a sin- 
gularity at 63U = 0, and the logarithm generates a branch cut along | arg(03?i)| = vr. The 
It-independent term 7e is the Euler-Mascheroni constant 7e = e~^ In x dx and is obtained 
by expanding e~'^^ in series, integrating term-by-term, and matching the undetermined con- 
stant of integration with the left-hand side of (B.4) in the limit n — )• 0. 

Because (B.4) is uniformly convergent we may exchange integration and summation, 
evaluating Iq using term-by-term integration. We find [34] 



loiks) 



n=l 



(B.5) 



This is singular when i?3 — )• 0, which corresponds to the squeezed limit where either ki or 
k2 — >■ 0. There is nothing unphysical about this arrangement of momenta, for which 63 — t- 00, 
but we will encounter a power-law divergence unless the bracket [• • • ] vanishes sufficiently 
rapidly in the same limit. 

The sum converges absolutely if —1/2 < Re(^3) < 1/2, corresponding to the narrow 
physical region < k^ < kt/4: where ^3 is being squeezed. For 9^ satisfying this condition 
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the summation can be performed explicitly, yielding ln(l + 2^3)"^ = lni?3. It can be verified 
that 10(^3) is analytic in the half-plane Re(03) > 0. Therefore it is possible to analytically 
continue to 63 ^ 1/2, allowing the behaviour as ^3 — >• 00 and '!?3 — )• to be studied. We find 

loih) = -— ^ ln(l - ^3) = --^Jo{k3), (B.6) 
VsktCs ktCs 

where we have introduced a function Jq, which is related to Iq by a simple change of normal- 
ization. It is more convenient to express our final bispectra in terms of Jq rather than Iq. In 
the limit 'd^ — )• 0, Eq. (B.6) is finite with Jq — )• —1. The ^3 — )• 00 behaviour of the sum has 
subtracted both the power-law divergence of the prefactor and the logarithmic divergence of 
ln(203). To obtain this conclusion required information about the entire momentum depen- 
dence of (B.l) in order to perform a meaningful analytic continuation.^^ Had we truncated 
the ksC dependence of (B.l), or equivalently the O^u dependence of (B.4), we would have 
encountered a spurious divergence in the squeezed limit and a misleading prediction of a 
strong signal in the local mode. Several of our predictions, including recovery of the Malda- 
cena limit described by (5.3), depend on the precise numerical value of Jq as 6*3 — 00 and 
therefore constitute tests of this procedure. 

We can now evaluate Imiks) for arbitrary m, although only the cases m ^ 2 are required 
for the calculation presented in the main text. The necessary expressions are 

In the limit ??3 — )• one can verify that Ji — >• —1/2 and J2 — )■ —2/3. 



C Useful integrals 



To simplify evaluation of the various integrals which arise in computing the bispectrum, it is 
helpful to have available a formula for a master integral, J^, for which the integrals of interest 
are special cases [34]. We define 



J. 



J — ( 



de e 



(C.l) 



^''A precisely analogous situation occurs at finite temperature, where one first computes the imaginary 
time Green's function G{iijJn) on a discrete spectrum of Matsubara frequencies {ujn}- To obtain the real- 
time Green's function one must extend G into the complex plane (in the sense that one finds a function 
G{z), where z £ C, which coincides with the imaginary time Green's function when z = icj„), and analytically 
continue to the real axis. Unless one has knowledge of G on the complete set of Matsubara frequencies — which 
may be impossible in an effective field theory — it is difflcult to extract meaningful results for the real-space 
correlation function. In the present case, the requirement to know the series expansion in (B.4) or (B.5) for 
all n is analogous to the requirement to know G(ia;„) for all Matsubara frequencies. 
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where A^^, represents the number of e- folds elapsed since horizon exit of the reference mode k^,, 
giving = In \k-^Csi,^\. Eq. (C.l) is an oscillatory integral, for which asymptotic techniques 
are well-developed, and which can be evaluated using repeated integration by parts. To 
evaluate the A^*-dependent terms one requires the standard integral 



lim 



d^ A^*e 



ikfCs 



he 



1 / vr 
7E + iTT 



(C.2) 



which can be obtained by contour rotation and use of the Euler-Mascheroni constant defined 
below Eq. (B.4). We conclude 

1 



hcs 



70 



71 + 5i 272 + 3(52 , 673 + 11^3 _^ 2474 + 50 J4 



kf 



kl 



kf 



We use this result repeatedly in §3.3 to evaluate integrals in closed form. 



(C.3) 



D The special functions R and Q 

For completeness, we give explicit formulae for the functions Q{ki, k2, ks) and R{ki, k2, k^) 
which were used to express a subset of the slow-roll corrections in Ref . [43] . These functions 
can be written in closed form using the integrals Jq, Ji and J2 defined in Appendix B. To 
do so, we introduce a quantity a(ki) which satisfies 

kt 



a{ki) = 



1. 



(D.l) 



We are adopting our usual convention in which a{ki) is a function of all ki, but only the 
asymmetrically occurring momentum is written explicitly. Also, in this Appendix, the Greek 
symbols a, f3 and 7 are used independently of their meaning elsewhere in this paper. The 
combination can be written [cf. (B.2)] 

a{ki) — 1 



' aiki) + 1 

The definitions of R and Q can be found in Eqs. (B.84)-(B.86) of Ref. [43]. We find 



R 



'^2"'3 

"fcT 



[l + a(A;i)]2 l + a{ki) 



Mki) - Ji{ki) + 



a{ki)J2{ki 
[l + a{ki)] 



+ sym. , 



(D.2) 



(D.3) 



where 'sym.' denotes the symmetric sum formed by adding the two distinct combinations 
kl — >• k2 and /ci — ?■ ^3. Eq. (D.3) is accurate up to terms which vanish like powers of |fcr|, and 
are therefore negligible a few e-folds after horizon crossing. With the same understanding, 
we find that Q can be written 



Q 



Hki 



(4 + 2[a{ki) - 1][7E + In 2]) + {kl + 

1 



5o(A;i)(ln[l + a(A;i)]+7E-l) + 



l + a{ki] 



In — + 7E 
kl 



h{ki)Uki) 



72(^1 



1 + a[ki 
Mki)Ji{ki) 



l + a{ki) 



+ sym. 
(D.4) 
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The functions appearing here are 



/5o(fci) = -|E^' (D.5) 



6 

i 



i 

7i(A:i) = -h{kl + kl) - kik2ksa{ki) - ^^^^ ^ fcf ^ 

i 

72(fei) = -A;ifc2fc3a(fci) - - V (D.IO) 

5o(fci) = [7o(A;i)+/3o(A:i)][l + a(A:i)] + 7i(A;i). (D.ll) 
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